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ABSTRACT 


This  thesis  presents  the  design  of  a  model  reference  adap¬ 
tive  control  system  for  a  double  effect  evaporator.  Liapunov's 
second  or  direct  method  of  stability  analysis  is  used  in  the  design 
approach  to  ensure  a  stable  system  response.  Hybrid  computer  simula¬ 
tion  runs  were  performed  in  order  to  test  the  effect  of  several  design 
parameters  on  the  control  scheme.  Experimental  runs  were  also  con¬ 
ducted  using  an  IBM  1800  process  control  computer  and  the  pilot  plant 
process . 

A  literature  survey  of  hybrid  computer  applications  in  the 
chemical  engineering  field  is  included.  An  EAI  590  hybrid  computer 
system  was  recently  installed  in  the  Department  of  Chemical  and 
Petroleum  Engineering  and  the  literature  survey  reveals  several  pos¬ 
sible  applications  for  this  hybrid  computer  system.  The  literature 
survey  is  accompanied  by  several  example  problems  that  demonstrate 
hybrid  computer  techniques. 

The  existing  theory  for  model  reference  adaptive  control 
designed  by  Liapunov's  second  method  is  presented  for  a  linear  time- 
invariant  process  with  a  basic  feedforward,  feedback  control  config¬ 
uration.  Extensions  to  the  basic  theory  for  adaptive  integral  and 
setpoint  control  are  also  included.  The  basic  feedback  control 
design  is  then  applied  to  a  fifth-order  linear  state-space  model  of 
a  double  effect  evaporator  using  an  optimal  feedback  reference  model. 

Hybrid  simulation  results  demonstrate  the  ability  of  model 
reference  adaptive  control  to  improve  upon  the  initial  control  policy 
and  eventually  attain  the  same  response  as  the  reference  model.  The 
control  scheme  is  able  to  survive  unmeasured  disturbances,  poor 
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initial  control  policies,  and  physically  impossible  reference  models. 
The  effect  of  other  design  parameters  on  control  is  also  demonstrated. 
The  control  scheme  demonstrates  the  ability  to  adjust  and  improve 
control  with  continuing  process  disturbances.  However,  numerical 
problems,  associated  with  the  discrete  approximation  to  the  continuous 
adaptive  control  law,  develop  with  large  adaptive  loop  gains  or  large 
control  intervals.  Low  adaptive  loop  gains  and  a  well-tuned  multiloop 
initial  control  policy  produce  excellent  control  and  model  following 
for  major  process  disturbances. 

The  experimental  runs  that  test  the  model  reference  adaptive 
control  system  verify  the  simulation  results  and  demonstrate  the 
industrial  feasibility  of  this  control  scheme. 
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CHAPTER  ONE 
INTRODUCTION 

Automatic  control  systems  inevitably  require  on-line  tuning 
to  achieve  a  satisfactory  degree  of  control.  Furthermore,  frequent 
readjustment  of  controller  constants  may  be  required  if  process  condi¬ 
tions  change  significantly.  An  adaptive  control  system  is  one  in  which 
the  control  parameters  or  control  configuration  changes  automatically 
during  process  operation  to  compensate  for  changing  process  conditions. 
In  the  model  reference  approach  to  adaptive  control,  the  control  para¬ 
meters  are  adapted  so  that  the  process  response  approaches  that  of  a 
reference  model  when  both  are  subjected  to  the  same  inputs. 

Modern  multivariable  control  systems  often  rely  heavily  on 
model  accuracy  for  their  success.  However,  in  many  industrial  situa¬ 
tions  a  process  model  may  not  be  readily  available.  Two  adaptive 
control  approaches  reduce  the  requirements  for  model  accuracy:  model 
reference  adaptive  control  (MRAC)  and  optimal  adaptive  control.  The 
reference  model  in  MRAC  represents  the  desired  closed-loop  plant 
dynamics  rather  than  the  actual  plant  dynamics  and  thus  does  not  require 
the  availability  of  an  accurate  process  model.  Optimal  adaptive  control 
combines  process  identification  with  optimization  of  the  continually 
updated  process  model.  MRAC  has  the  advantage  of  not  requiring  process 
identification  and  time  consuming  optimization  calculations  and  would 
seem  to  be  more  applicable  to  industrial  situations.  Only  MRAC  is  dealt 
with  in  this  thesis. 

MRAC  schemes  have  found  extensive  use  in  the  aerospace  indus¬ 
tries  both  in  single-input,  single-output  and  multivariable  control 
systems.  However,  the  process  industries  have  only  applied  MRAC  to 
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single  variable  systems.  This  thesis  is  concerned  with  the  application 
of  multivariable  MRAC  to  a  pilot  plant  double  effect  evaporator. 

The  pilot  plant  double  effect  evaporator  in  the  Department  of 
Chemical  and  Petroleum  Engineering  at  the  University  of  Alberta  has  been 
the  subject  of  several  studies  in  modelling  and  control  [1  -  4].  The 
experimental  verification  of  many  conventional  and  advanced  control 
schemes  has  been  made  possible  by  interfacing  the  evaporator  to  an  IBM 
1800  process  control  computer. 

Early  studies  in  the  Department  have  investigated  multi  loop 
Direct  Digital  Control  (DDC)  including  inferential  and  feedforward 
control  configurations  [2,3].  Recent  work  has  dealt  with  optimal 
regulatory  multivariable  control  [4]  and  control  techniques  for  making 
optimal  setpoint  changes  [5].  These  control  techniques  proved  much 
more  successful  than  conventional  schemes.  Current  research  [6]  deals 
with  control  techniques  using  reduced  order  state-space  models. 

Studies  are  also  being  conducted  on  multivariable  control  techniques 
[7,8]  which  consider  various  conditions  (noise,  time-delays)  that  are 
not  significant  problems  for  the  evaporator  but  can  be  important  in 
other  processes.  However,  in  all  these  studies  an  approximate  dynamic 
model  of  the  evaporator  was  required. 

Work  on  a  nine  inch,  eight  tray  distillation  column  in  the 
Department  [9]  has  illustrated  the  difficulty  and  time  involved  in 
obtaining  an  adequate  state-space  model  for  some  systems.  Part  of  the 

t 

motivation  for  this  study  was  to  evaluate  an  adaptive  control  technique 
that  relies  less  heavily  on  an  accurate  model  of  the  process  being 
controlled. 

This  thesis  presents  a  study  of  a  MRAC  technique  designed  via 
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Liapunov's  direct  method  of  stability  analysis.  The  MRAC  technique, 
based  on  an  approach  developed  by  Winsor  and  Roy  [11]  and  Porter  and 
Tatnall  [10],  was  first  evaluated  by  using  a  hybrid  simulation  of  the 
IBM  1800/evaporator  configuration.  The  simulation  results  were  then 
experimentally  verified  on  the  pilot  plant  system.  Since  this  thesis 
research  was  the  first  project  in  the  Department  of  Chemical  and 
Petroleum  Engineering  to  use  the  recently  installed  EAI  590  hybrid 
computer  system,  several  hybrid  programs  demonstrating  various  aspects 
of  hybrid  computation  were  developed  in  order  to  become  familiar  with 
the  system. 

1.1  OBJECTIVES  OF  THE  STUDY 

The  intent  of  this  thesis  is  to  experimentally  verify  a  model 
reference  adaptive  control  technique  designed  via  Liapunov's  direct 
method.  The  theory  for  this  adaptive  control  design  approach  was 
developed  in  the  late  1960's  [10,11]  and  has  been  investigated  in 
several  mathematical  and  simulation  studies.  However,  the  only  experi¬ 
mental  verification  of  the  simulation  results  has  been  conducted  on  a 
simple  hydromechanical  process  modelled  as  a  first-order  system  [12]. 

The  following  steps  were  taken  to  achieve  the  completion  of 
this  project: 

i)  Develop  familiarity  with  the  new  hybrid  computing  system 
through  the  solution  of  typical  problems  taken  from  the 
literature. 

ii)  Design  a  MRAC  system  for  the  evaporator  using  Liapunov's 
direct  method. 

iii)  Simulate  the  evaporator  control  system  using  the  EAI  590 
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hybrid  system  and  use  the  hybrid  simulation  to  investi¬ 
gate  the  critical  design  features  of  this  method, 
iv)  Experimentally  verify  the  simulation  results  using  the 

actual  double  effect  evaporator  and  the  IBM  1800  computer. 

1.2  STRUCTURE  OF  THE  THESIS 

The  thesis  involves  two  distinct  sections:  the  first  section 
presented  in  the  next  chapter  includes  a  literature  survey  of  representa¬ 
tive  hybrid  computer  applications  in  chemical  engineering  and  illustrates 
some  of  the  applications  on  the  EAI  590  system  in  a  series  of  sample 
problems.  The  problems  were  used  to  test  the  capabilities  of  the  analog, 
digital,  and  interface  sections  of  the  hybrid  system.  They  also  pro¬ 
vided  insight  into  which  programming  language  would  be  best  for  the 
model  reference  adaptive  control  simulation  to  follow. 

The  second  section  of  the  thesis  involves  MRAC.  Chapter  Three 
presents  a  literature  survey  of  MRAC  techniques  that  have  developed 
over  the  last  ten  to  fifteen  years.  The  theory  of  the  MRAC  design 
approach  via  Liapunov's  direct  method  and  the  specific  application  to 
the  double  effect  evaporator  are  then  presented. 

The  hybrid  simulation  results  are  given  in  Chapter  Four  and 
conclusions  regarding  the  usefulness  of  MRAC  are  drawn  from  these 
results.  Chapter  Five  presents  the  experimental  verification  of  the 
simulation  results  for  the  actual  evaporator  using  the  IBM  1800  process 
control  computer.  The  experimental  runs  also  produced  some  unexpected 
problems  not  encountered  in  the  simulation  studies.  The  computer  pro¬ 
grams  used  to  obtain  the  simulation  and  experimental  results  are  given 
in  detail  in  research  reports  and  user's  manuals  that  are  referenced  in 
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each  chapter. 

Chapter  Six  summarizes  the  overall  conclusions  drawn  from  the 

study. 
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CHAPTER  TWO 

HYBRID  COMPUTER  APPLICATIONS  IN  CHEMICAL  ENGINEERING 
2.1  INTRODUCTION 

Hybrid  computers  have  found  increasing  use  in  chemical  engineer¬ 
ing  simulations  since  the  early  1960's.  Many  problems  in  the  fields  of 
process  dynamics,  process  control,  and  optimization  are  best  suited  to 
a  hybrid  computer  because  a  hybrid  system  combines  the  high  speed 
integration  capabilities  of  an  analog  computer  with  the  arithmetic  and 
logic  capabilities  of  a  digital  machine.  The  improvements  in  digital 
software  and  interface  hardware  in  the  last  ten  years  have  helped  to 
increase  the  popularity  of  hybrid  computation  and  have  led  to  highly 
complex  simulations  of  large  systems  [1]. 

With  the  recent  purchase  of  an  EAI  590  hybrid  computer  system 
by  the  Department  of  Chemical  and  Petroleum  Engineering  of  the  University 
of  Alberta,  an  investigation  was  undertaken  in  order  to  determine  the 
capabilities  and  possible  applications  of  the  system.  This  investiga¬ 
tion  took  the  form  of  a  literature  survey  and  a  series  of  example  pro¬ 
blems  solved  on  the  EAI  590  system.  The  objectives  of  this  investiga¬ 
tion  were  as  follows: 

i)  To  determine  the  type  of  problems  applicable  to  hybrid 
solutions  in  the  fields  of  chemical  engineering  and 
process  control . 

ii)  To  test  the  level  of  problem  complexity  that  could  be 
solved  on  the  relatively  small  EAI  590  system. 

iii)  To  investigate  the  types  of  software  languages  avail¬ 
able  on  the  EAI  590  and  determine  which  language  would 
be  most  useful  in  later  research. 
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iv)  To  provide  familiarity  with  the  new  computer  by 
"debugging"  sample  problems. 

This  chapter  briefly  reviews  the  structure  of  a  hybrid  com¬ 
puter  system  followed  by  a  literature  survey  of  hybrid  computer  appli¬ 
cations  in  chemical  engineering.  Finally,  seven  hybrid  sample  problems 
are  presented  in  order  to  illustrate  some  of  the  applications  discovered 
in  the  literature  survey. 

2.2  HYBRID  COMPUTER  STRUCTURE 

Figure  2.1  presents  a  block  diagram  of  a  typical 
hybrid  computing  system. 


ANALOG  INTERFACE  DIGITAL 


FIGURE  2.1  HYBRID  COMPUTING  SYSTEM  BLOCK  DIAGRAM 

The  specifications  for  the  present  EAI  590  configuration  at  the  Univer¬ 
sity  of  Alberta  are  presented  in  Table  2.1.  As  the  table  indicates, 
the  hybrid  computer  includes  an  EAI  640  digital  computer  with  16  K  of 
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TABLE  2.1 

EAI  590  HYBRID  COMPUTING  SYSTEM 
(October  1971) 


ANALOG  (580) 

1.  36  High  gain  amplifiers 

(12  integrators). 


2.  Repetitive  operation 
capabilities  to  .2 
Msec. 


3.  33  Potentiometers 
(25  digitally  set) . 

4.  1  Analog  logic  panel 
with  16  AND  gates,  2 
registers,  1  counter, 
2  timers,  4  compari- 
tors,  2  function 
switches . 


For  more  details  see 
E.A.I.  590  ANALOG/HYBRID 
Computing  System,  Pub. 

No.  008002055.0  March  1968. 


INTERFACE  (693) 

1 .  12  Analog  to  Digital 

Converters  (ADC) 
13  bit  conversion 
(33  u  sec. ) . 

2.  4  Digital  to  Analog 
Multipliers  (DAM) 
(zero  order  hold)  14 
bit  accuracy,  9.0  y 
sec.  conversion. 


3.  4  sense  lines, 

4  control  lines. 

4.  2  interrupt  lines 
and  an  ADC  finished 
interrupt. 


See  E.A.I.  693  II 
INTERFACE  SYSTEM,  Pub. 
No.  008003052-0, 
October  1970. 


DIGITAL (640) 

1.  16K,  1.65  y 
sec.  core  cycle 
time. 


2.  Teletype,  cas¬ 
sette  and  high 
speed  paper 
tape  input- 
output  systems. 
Expansion  to 
line  printer 
and  card  reader 
possible. 

3.  12  hardware, 

16  bit  registers. 

4.  64  instructions 
and  64  external 
interrupts . 


5.  Software  in¬ 
cludes  FORTRAN 
IV  compiler, 
ASSEMBLER  and 
HOI  ,.  Hytran 
Operations 
Interpreter. 

See  E.A.I.  640 
Reference  Handbook, 
Pub.  No.  8009118-0, 
June  1969. 
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core  storage  and  without  high  speed  bulk  storage  devices  (e.g.  disks)  or 
high  speed  input-output  devices.  Because  of  the  small  core  storage  and 
the  slow  peripherals,  the  EAI  640  can  be  classified  as  a  "mini -computer". 
The  EAI  580  analog  computer  is  relatively  small  and  is  limited  by  a 
shortage  of  operational  amplifiers.  The  functions  usually  associated 
with  each  of  the  hybrid  sections  will  now  be  outlined, 

i)  Analog 

The  basic  function  of  the  analog  computer  in  hybrid  simulations 
is  to  solve  the  dynamic  equations  that  arise  in  modelling  a  chemical 
process  by  integrating  with  electronic  analog  components.  Other  opera¬ 
tions  in  addition  to  integration  may  be  performed  such  as  summation, 
multiplication,  function  generation,  and  limited  logic  operations; 
however,  the  main  advantage  of  analog  simulation  is  high  speed,  accurate, 
continuous  integration. 

ii)  Digital 

The  digital  computer  performs  the  arithmetic  and  logic  func¬ 
tions  in  hybrid  simulations.  It  also  has  far  greater  storage  and  data 
playback  capabilities  than  does  the  analog.  Most  of  the  input-output 
is  done  through  the  digital  peripherals. 

iii)  Interface 

The  interface  controls  the  communications  between  the  discrete 
digital  section  and  the  continuous  analog  section.  Analog  to  digital 
converters  (ADCs)  allow  the  digital  computer  to  monitor  analog  voltages. 
Digital  to  analog  converters  ( DACs )  transform  a  digital  number  into  a 
voltage.  Also,  logic  control  lines  provide  communication  of  logic  sig¬ 
nals  between  the  analog  and  digital  sections  allowing  interrupt  handling, 
timing  control,  mode  control,  contact  closures,  potentiometer  settings, 
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etc.  The  interface  conversion  times  are  often  crucial  to  successful 
hybrid  simulations  and  rapid,  accurate  components  are  needed. 

2.3  LITERATURE  SURVEY 

Hybrid  computer  solutions  can  be  classified  according  to  the 
requirements  placed  on  the  analog  and  digital  sections  during  the  solu¬ 
tion.  The  category  —  analog-oriented,  balanced,  or  digital-oriented  -- 
a  problem  falls  into  is  governed  by  the  type  of  problem  and  the  capa¬ 
bilities  of  the  specific  hybrid  system.  Analog-oriented  problems  can 
use  the  digital  computer  to  check  and  initiate  a  pure  analog  simulation 
and  to  monitor  the  analog  voltages  for  later  analysis.  Balanced  pro¬ 
blems  usually  utilize  the  analog's  integration  capabilities  and  the 
digital's  arithmetic  capacity  to  generate  functions  not  easily  pro¬ 
grammed  on  the  analog  computer.  Digital-oriented  problems  use  the 
digital  computer's  ability  for  mode  and  timing  control  to  utilize  the 
analog  essentially  as  a  subroutine  to  solve  dynamic  equations. 

The  following  hybrid  literature  survey  begins  with  analog- 
oriented  problems,  demonstrates  many  balanced  applications  and  con¬ 
cludes  with  digitally-oriented  problems.  The  emphasis  is  on  chemical 
engineering  applications,  although  many  similar  examples  of  hybrid 
simulations  in  other  fields  (bio-medical,  aerospace,  communications, 
nuclear,  etc.)  can  be  readily  discovered  in  the  literature. 

2.3.1  Analog-oriented  Problems 

The  first  stage  beyond  pure  analog  computation  in  hybrid 
simulations  is  the  digital  set-up  and  check  of  analog  circuits. 

Greenberg  [2]  uses  computer  checkout  and  analog  debugging  as  the  first 
part  of  a  four  part  program  to  teach  hybrid  computation  to  undergraduates. 
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Babushkin,  Kogan  and  Ftybashov  [3]  indicate  the  digital  program  used  to 
aid  in  analog  set-up  may  be  able  to  perform  several  functions.  They 
outline  programs  that  are  able  to  reduce  the  dynamic  equations  to  a 
series  of  first-order  differential  equations,  provide  the  optimal 
scale  factors  for  the  analog  simulation,  and  automatically  set  the 
potentiometers  from  these  results. 

The  checkout  facility  usually  includes  a  static  check,  i.e. 
checking  the  analog  patchboard  by  placing  voltages  at  various  loca¬ 
tions  in  the  patch  panel  and  checking  component  outputs.  The  digital 
computer  can  also  perform  a  dynamic  check  of  the  analog  solution  by 
solving  the  differential  equations  numerically.  Hambury  [4]  states 
that  simple  interpretive  languages  are  ideal  for  this  type  of  hybrid 
application  since  the  ability  to  easily  change  the  digital  program  is 
more  important  than  the  need  for  rapid  digital  execution.  The  first 
illustrative  example  in  section  2.4  demonstrates  the  use  of  a  low  level 
interpretive  language,  the  Hytran  Operations  Interpreter  (HOI)  [5], in 
checking  the  analog  simulation  of  the  Sutro  Weir  experiment  used  in 
the  undergraduate  laboratory  at  the  University  of  Alberta. 

2.3.2  Digitally  Generated  Functions 

Hybrid  computers  have  found  extensive  use  in  solving  nonlinear 
ordinary  differential  equations  that  often  arise  when  modelling  the 
dynamics  of  a  chemical  process.  Analog  nonlinear  elements  are  expensive 
and  inaccurate,  and  because  of  this,  many  analog  simulations  involve 
linearized  models.  Digital  computers  have  the  arithmetic  capabilities 
to  generate  highly  nonlinear  functions;  thus  digital  simulations  for 
nonlinear  problems  are  popular.  The  digital  computer  in  a  hybrid 
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simulation  performs  the  nonlinear  function  generation  and  the  analog 
integrators  eliminate  some  of  the  problems  (round-off  error,  trunca¬ 
tion  error,  and  computation  time)  that  can  develop  with  pure  digital 
simulations.  However,  the  synchronization  of  the  analog  and  digital 
computers  and  conversion  errors  (analog  to  digital  and  digital  to 
analog  conversions)  add  extra  problems  to  hybrid  simulations. 

The  major  drawbacks  to  digitally  generated  functions  in 
hybrid  simulations  arise  from  sampling  a  continuous  analog  function. 

The  analog  is  sampled  at  a  time,  t,  and  the  function  generated  at  a 
time,  t  +  At,  where  At  is  the  time  required  to  calculate  the  function 
and  complete  the  conversions  required.  This  causes  the  function  to  be 
in  error  by  the  amount  of  time  required  to  complete  the  function 
generation.  Rapid  ADCs  and  DACs  can  make  this  error  negligible  in 
most  real-time  and  moderately  time-scaled  (accelerated)  simulations. 

When  the  analog  is  operating  in  the  repetitive  operation  mode  or 
simulating  a  rapidly  changing  process,  these  errors  can  become  critical. 
The  text  by  Bekey  and  Karplus  [6]  presents  a  complete  chapter  (Chapter 
5)  on  errors  caused  by  these  conversions.  Methods  of  predicting  the 
value  of  a  function  at  time,  t  +  At,  from  an  input  at  time,  t,  have 
been  extensively  investigated  [6,7].  Vansteenkiste  [8]  has  summarized 
some  of  these  compensating  formulae  for  both  zero-order  and  first-order 
hold  circuits  that  are  present  in  the  DACs  of  most  hybrid  systems. 
Methods  of  testing  these  formulae  have  also  been  developed  [9,10].  The 
hybrid  example  problem  in  section  2.4  demonstrates  the  accuracy  of 
digitally  generated  functions  by  comparing  a  digitally  generated  sine 
wave  to  an  analog  sine  wave  for  different  sampling  frequencies. 

Many  examples  of  digitally  generated  functions  can  be  found  in 
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hybrid  simulations  of  process  systems.  Severe  nonlinearities  occur  in 
simulations  of  chemical  reactor  systems  where  the  reaction  rate  involves 
the  Arrhenius  relationship,  expC-E/RT),  and  is  not  suitable  for  linear¬ 
ization.  Sharpe  [11],  Dahlin  and  Nelson  [12],  and  Ruszkay  and  Mitchell 
[13]  all  demonstrated  the  use  of  a  digital  computer  in  nonlinear  func¬ 
tion  generation  for  hybrid  solution  of  reaction  problems.  An  example 
of  the  hybrid  simulation  of  a  nonlinear  batch  reactor  is  presented  in 
the  example  problem  section.  Sharpe's  article  [11]  contains  a  compari¬ 
son  of  the  time  and  cost  for  a  hybrid  simulation  and  a  digital  simula¬ 
tion  of  space  vehicle  dynamics.  A  similar  comparison  is  made  in  the 
reactor  problem  presented  in  the  next  section.  Both  simulations 
incurred  large  savings  on  the  hybrid  computer. 

A  mathematical  function  ideally  suited  to  digital  generation 
is  a  time  delay.  Pure  analog  simulation  of  time  delays  involves  a 
series  approximation  (Pade,  Stubbs-Single  etc.)  that  requires  many 
operational  amplifiers  and  often  yields  good  results  for  only  certain 
types  of  simulations.  By  contrast,  simulating  time  delays  on  a  digital 
machine  is  very  accurate.  Fox,  Parsons  and  Dutcher  [14]  simulate  time 
delays  digitally  in  the  hybrid  simulation  of  a  polymerization  process. 
Nil  son  [15]  has  demonstrated  a  technique  for  generating  delayed  func¬ 
tions,  f(t-i),  from  a  variety  of  functions,  f(t).  Section  2.4  presents 
examples  of  time  delay  programs. 

2.3.3  Sampled-Data  Systems 

The  sampling  problems  inherent  in  hybrid  simulations  of  con¬ 
tinuous  systems  can  be  put  to  advantage  in  analysing  sampled-data 
systems.  As  Figure  2.2  indicates,  the  analysis  of  sampled-data  systems 
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allows  for  a  sample  and  hold  device  present  in  the  system. 


FIGURE  2.2  SAMPLED-DATA  CONTROL  SYSTEM 

For  sampled-data  problems  analysed  in  the  Z-domain  [16],  the 
digital  computer  acts  as  the  digital  controller,  the  interface  com¬ 
ponents  represent  the  sample  and  hold  device,  and  the  analog  computer 
simulates  the  continuous  process.  A  digital  controller  designed  to 
give  a  dead-beat  response  [17]  is  demonstrated  in  section  2.4. 

Digital  control  of  a  process  always  requires  sampling.  Often 
the  effect  of  sampling  is  ignored  and  the  controlled  process  is  analysed 
as  a  continuous  system.  Many  digital  simulations  do  not  allow  for  errors 
that  can  arise  because  of  sampling.  However,  if  the  analog  simulates 
the  actual  process  and  the  digital  controls  this  process,  sampling 
errors  may  be  discovered  at  the  time  of  the  simulation,  not  when  the 
control  scheme  is  being  implemented.  Fox  et  al  [14]  used  a  hybrid  com¬ 
puter  to  simulate  Direct  Digital  Control  (DDC)  of  a  polymerization 
process  and  investigated  the  effect  of  sampling  frequency.  Lawrence  et 
al  [18]  used  the  similarity  between  hybrid  simulation  and  control  of 
an  actual  plant  to  investigate  DDC  control  algorithms  on  control  loops 
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that  existed  at  a  nuclear  power  station.  More  examples  of  hybrid 
digital  control  simulations  will  be  reviewed  in  the  next  section  on 
control  applications. 

2.3.4  Control 

The  ability  of  the  analog  to  provide  rapid  solutions  to 
dynamic  equations  has  caused  extensive  use  of  hybrid  computers  in 
control  studies.  Single  variable  control  problems  have  been  studied 
on  hybrid  computers  to  obtain  optimal  control  settings  for  conventional 
controllers.  Sharpe  [11]  describes  a  hybrid  simulation  of  a  tubular 
chemical  reactor  used  to  optimize  the  controller  constants  for  a  PID 
controller.  Ung  [20]  designed  a  three  mode  controller  for  a  heat 
exchanger  control  system  by  minimizing  the  integral  of  the  absolute 
error  performance  criterion.  A  similar  optimization  problem  for  the 
Sutro  Weir  experiment  is  presented  in  section  2.4. 

Multivariable  control  studies  have  been  popular  hybrid  com¬ 
puter  applications  because  the  implementation  of  multivariable  schemes 
depends  on  a  simultaneous  knowledge  of  all  the  controlled  variables 
(available  because  of  the  analog  parallel  integration  capabilities)  and 
on  matrix  operations  (available  from  the  arithmetic  capabilities  of 
the  digital).  Bekey  [21]  has  used  hybrid  computers  in  the  dynamic 
optimization  of  multivariable  systems.  As  Figure  2.3  illustrates, 
the  analog  simulates  the  process  represented  by  vector-matrix  differen¬ 
tial  equations.  The  digital  evaluates  the  criterion  to  be  minimized 
from  the  analog  results  (controlled  variables,  2l)  and  calculates  the 
control  action,  u,  from  a  control  algorithm  depending  on  the  criterion, 
the  controlled  variables,  and  the  control  action. 
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FIGURE  2.3  HYBRID  COMPUTERS  IN  MULTIVARIABLE  CONTROL  SYSTEM  DESIGN 

Gonzalez  [22]  has  used  the  same  hybrid  approach  to  calculate  the  mini¬ 
mum  time  optimal  control  of  hypothetical  second-and  third-order  pro¬ 
cesses  with  constrained  control  variables. 

Hybrid  computers  have  found  extensive  use  in  the  study  of 
adaptive  control  techniques.  Sandoz  and  Swanick  [23]  conducted  a 
hybrid  simulation  of  a  model  reference  adaptive  control  system  in  which 
the  error  between  the  analog  "process"  and  digitally  generated  model 
was  minimized  by  adapting  the  control  of  the  system  at  each  interval. 
They  reported  excellent  results  for  a  second-order  system  with  feedback 
control.  Godebole  and  Smith  [24]  tested  another  adaptive  control  tech¬ 
nique  using  both  digital  and  hybrid  simulations.  The  control  scheme 
greatly  improved  control  for  systems  modelled  with  transfer  functions. 
Furthermore,  they  maintain  the  hybrid  simulation  to  be  the  most 
realistic  because  the  analog  process  simulation  and  the  digital  control 
calculations  are  working  in  parallel;  therefore  sampling  errors  that 
would  occur  in  actual  process  control  also  occur  in  the  simulation. 

By  contrast,  digital  simulation  suspends  process  dynamics  until  the 
control  calculations  are  complete.  The  later  chapters  of  this  thesis 
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illustrate  a  model  reference  adaptive  control  technique  simulated  on 
the  EAI  590  hybrid  system. 

2.3.5  Real-Time  Applications 

Real-time  applications  of  hybrid  computers  are  a  relatively 
recent  development  but  could  have  important  applications  in  the  future. 
Mel li champ  [25,26]  demonstrated  this  application  when  he  optimized  a 
second-order  process  with  [25]  and  without  [26]  time  delays.  He 
optimized  a  process  model  on  the  analog  in  "fast-time"  and  returned  the 
optimal  control  values  to  the  process  in  real-time.  Then,  the  process 
states  were  read  by  the  digital  control  computer  and  used  to  update  the 
analog  model  for  the  next  optimization  iteration.  Mellichamp  reported 
excellent  results.  With  process  control  computers  becoming  more  common¬ 
place,  this  type  of  hybrid  real-time  application  should  become  more 
popular. 

2.3.6  Optimization  and  Parameter  Estimation 

Optimization  and  parameter  estimation  problems  rely  heavily 
on  the  rapid  solution  of  an  analog  model.  Andrews  and  Moore  [19]  used 
a  hybrid  simulation  to  optimize  the  cost  of  the  catalytic  dehydro¬ 
genation  of  normal  butane.  The  analog  computer  simulated  the  reactor 
section  of  the  plant  and  contained  an  optimization  routine  to  minimize 
the  reactor  operating  cost.  The  digital  computer  updated  the  optimiza¬ 
tion  parameters  in  the  analog  section,  calculated  the  distillation 
costs  from  the  reactor  outlet  composition,  and  from  the  combined  total 
cost  (reactor  and  distillation)  optimized  the  entire  process  by  varying 
reactor  conditions  (feed  rates,  temperatures,  recycle  rate,  etc.). 

In  parameter  estimation,  the  digital  stores  the  actual  process 
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response  and  the  analog  produces  a  model  response.  The  parameters  of 
the  model  are  then  updated  until  the  desired  "goodness  of  fit"  between 
the  process  and  model  are  attained.  Carlson  [27]  outlines  the  pro¬ 
cedure  for  parameter  estimation  in  kinetic  reaction  equations.  His 
paper  is  typical  in  that  it  uses  the  rapid  analog  solution  of  the 
model  equations  to  quickly  obtain  the  best  parameters  (rate  constants) 
for  the  kinetic  rate  equations. 

2.3.7  Simulation  of  Staged  Processes 

The  simulation  of  staged  processes  is  a  major  digital - 
oriented  hybrid  application.  In  the  hybrid  simulation  of  staged  pro¬ 
cesses,  the  process  being  simulated  is  divided  into  identical  stages 
and  the  differential  equations  describing  the  dynamics  of  the  stage 
are  programmed  on  the  analog.  The  digital  then  links  the  simulation 
together  by  transferring  the  boundary  conditions  between  stages.  The 
stages  may  be  real,  as  in  a  plate  distillation  column,  or  contrived, 
as  in  a  packed  column  or  tubular  reactor.  For  some  processes,  a  large 
number  of  stages  are  required  for  an  accurate  model.  Therefore,  an 
excessive  number  of  components  would  be  required  in  pure  analog  simula¬ 
tion.  Furthermore,  accuracy  and  computation  time  problems  could  arise 
in  digital  simulations.  Many  examples  of  large  scale,  staged  process, 
hybrid  simulations  are  available  in  the  literature. 

Dahl  in  and  Nelson  [12]  used  a  hybrid  computer  to  study  the 
cooling  of  a  tubular  chemical  reactor  by  countercurrent  flow.  They 
modelled  the  tubular  reactor  as  a  series  of  backmix  reactors  with 
backflow  between  adjacent  stages.  Two  first-order  consecutive  reactions 
were  considered  involving  three  dynamic  component  balances  and  a  heat 
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balance  for  each  stage.  A  single  stage  analog  circuit  used  12  opera¬ 
tional  amplifiers,  and  if  20  stages  were  needed  to  model  the  reactor, 
then  over  200  operational  amplifiers  would  be  required.  Few  analog 
computers  have  this  capacity.  With  the  hybrid  solution,  each  stage 
was  solved  separately  and  the  solution  of  that  stage  was  read  into  the 
digital  via  ADC  channels.  These  values  were  then  used  to  calculate 
the  boundary  conditions  for  the  next  stage.  Therefore,  all  the 
identical  stages  used  the  same  analog  solution.  The  countercurrent 
cooling  of  the  reactor  led  to  a  split-boundary  value  problem  and  an 
iterative  solution  was  required.  The  outlet  coolant  temperature  had 
to  be  assumed  in  order  to  solve  the  problem  and  the  resulting  solution 
for  the  inlet  cooling  temperature  was  compared  to  a  known  inlet  condi¬ 
tion.  The  simulation  had  to  be  repeated  if  the  temperatures  were  not 
within  reasonable  agreement.  Because  of  the  iterative  nature  of  the 
solution,  the  rapid  single  stage  analog  solution  made  the  hybrid 
calculation  time  much  less  than  the  digital  solution  time. 

Frank  and  Lapidus  [28]  used  this  technique  to  simulate  a 
seven  stage,  five  component  distillation  column  and  a  fixed  bed 
catalytic  reactor  modelled  as  eight  continuous  stirred  tank  reactors. 
Ruszkay  and  Mitchell  [13]  simulated  a  twenty-five  stage  reacting 
distillation  column  using  the  single  analog  stage  technique.  This 
method  of  solution  is  equivalent  to  solving  a  partial  differential 
equation  in  length  and  time  coordinates  by  the  discrete-space-discrete¬ 
time  (DSDT)  method.  Indeed,  the  hybrid  solution  of  partial  differential 
equations  is  an  extremely  popular  hybrid  application. 
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2.3.8  Partial  Differential  Equations 

Partial  differential  equations  (PDE)  often  arise  when  model¬ 
ling  chemical  processes.  When  analytic  solutions  are  not  available, 
the  PDE  can  be  solved  using  numerical  techniques  on  digital  computers, 
however,  large  amounts  of  computer  time  are  involved  in  multidimensional 
integrations.  In  the  past  few  years,  several  hybrid  techniques  for 
solving  PDE  have  been  investigated;  continuous-time-discrete-space(CTDS) , 
continuous-space-discrete-time  (CSDT),  and  discrete-space-discrete-time 
(DSDT).  These  methods  are  presented  in  detail  in  Bekey  and  Karp! us  [6] 
(Chapter  8).  Essentially  all  the  methods  involve  solving  one  dimension 
of  the  differential  equation  (time  or  space)  continuously  on  the  analog 
and  the  other  dimension  by  discrete  methods  on  the  digital.  The 
dynamics  of  many  systems  are  best  described  by  a  PDE  and  many  hybrid 
solution  techniques  have  been  discussed  in  the  literature. 

Carling  [29]  developed  equations  for  the  hybrid  simulation  of 
a  counterflow  heat  exchanger  with  thermal  capacity  in  both  fluids  and 
the  tube  wall.  He  was  unable  to  achieve  a  complete  hybrid  solution 
because  of  hybrid  hardware  difficulties  on  his  computer.  Carling  did 
verify  his  technique  on  a  digital  hybrid  simulator.  Poulsen  [30] 
solved  a  similar  problem  for  the  dynamics  of  steam  loss  from  a  tube. 

Owen  and  Kistler  [31]  and  Vichnevetsky  and  Tomalesky  [32]  have  developed 
two  different  techniques  to  obtain  the  heat  exchanger  solution.  The 
same  heat  exchanger  problem  described  by  Owen  and  Kistler  is  solved  in 
the  next  section  using  the  standard  CSDT  method. 

A  PDE  similar  to  those  arising  from  heat  exchanger  modelling 
is  found  in  fixed  bed  adsorption  models  and  has  been  solved  by  Eteson 
and  Zwiebel  [33]  on  a  hybrid  computer. 
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Second-order  partial  differential  equations  in  the  spatial 
coordinates  often  occur  when  diffusion  in  mass  transfer  or  conduction 
in  heat  transfer  is  important.  Bryant,  Amiot,  et  al  [42]  have  solved 
the  two  dimensional  heat  transfer  differential  equation  of  the  form 


1  _1  (x  ii) 

X  3X  v  3X; 


(2.1) 


where  x  is  radial  distance,  z  is  axial  distance  and  £  is  temperature. 
This  equation  results  from  a  steady-state  laminar  flow  heat  balance 
around  a  double  pipe  heat  exchanger.  Equation  (2.1)  is  a  form  of  the 
Sturm-Liouville  equation  and  often  occurs  with  split-boundary  values. 
The  application  of  numerical  techniques  to  the  hybrid  solution  of  this 
type  of  PDE  has  made  use  of  high  speed  analog  integration  to  reduce 
solution  time. 

The  diffusion  equation 


a  3t 


(2.2) 


has  been  solved  many  times  in  the  literature  [34-36].  Many  approaches 
to  solving  this  equation  have  been  discussed  because  at  high  solution 
speeds,  sampling  errors  can  be  significant.  Also,  difficulties  are 
encountered  in  finding  a  value  of  the  time  increment  (CSDT  method) 
that  yields  an  acceptable  compromise  between  the  effect  of  round-off 
error  (in  analog  to  digital  and  digital  to  analog  conversions)  and  the 
truncation  error  associated  with  time-discretization. 

One  of  the  most  complex  PDE  problems  solved  by  hybrid  tech¬ 
niques  to  date  deals  with  river  pollution  [37]  and  takes  the  form 
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If  -  £ k  £  -  ^  -  «*>  +  «*■*>  <2-3> 

where  v  =  river  velocity 

D(x)  =  pollutant  degradation 
C  =  pollutant  concentration 

k  =  diffusion  constant 

and  f(x,t)  =  pollutant  source  function 

Third-order  and  higher  partial  differential  equations  are  rarely  solved 
on  hybrid  computers,  mainly  because  they  seldom  occur  in  process 
modelling.  Multidimensional  (more  than  two  dimensions)  simulations 
also  do  not  occur  because  the  synchronization  of  two  analog  machines 
would  be  required  for  the  solution  unless  all  but  one  dimension  was 
solved  for  numerically. 

The  papers  just  reviewed  cover  all  the  major  areas  of  hybrid 
computer  problems  in  chemical  engineering  and  are  representative  of 
work  in  this  field.  From  the  literature  survey,  some  applications  of 
the  University  of  Alberta  EAI  590  hybrid  system  can  be  suggested.  These 
are: 

i)  Evaluation  of  Control  Techniques:  The  ability  to  adjust 
the  time  scale  and  the  close  correlation  of  hybrid  simulation  and  pro¬ 
cess  operating  conditions  under  computer  control  make  hybrid  computers 
ideal  for  testing  computer  control  techniques  before  actual  process 
implementation.  The  following  chapters  present  the  first  such  research 
performed  using  the  department's  EAI  590  hybrid  computer. 

ii)  Simulation  of  Staged  Processes:  The  solution  of  models 
of  present  equipment  at  the  University  of  Alberta,  such  as  the  8  tray. 
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9  inch  distillation  column  could  be  attempted  on  the  hybrid  computer. 

iii)  Solution  of  Partial  Differential  Equations:  The  partial 
differential  equations  describing  heat  exchangers  and  mass  transfer 
operations  could  be  solved  on  the  EAI  590  computer. 

iv)  Process  Identification  and  Kinetic  Parameter  Estimation 

v)  Undergraduate  Teaching  Aid:  Hybrid  computers  are 
becoming  more  commonplace  and  undergraduate  courses  in  hybrid  computer 
techniques  are  being  taught  at  some  universities  [2], 

2.4  ILLUSTRATIVE  EXAMPLES 

To  demonstrate  some  of  the  hybrid  applications  mentioned  in 
the  literature  survey,  several  sample  problems  were  solved  on  the  EAI 
590  hybrid  computer.  The  seven  example  problems  are  the  following: 

2.4.1  Digital  Initiation  and  Check  of  an  Analog  Simulation 

2.4.2  Reconstruction  of  Sampled  Signals 

2.4.3  Simulation  of  Time  Delays 

2.4.4  Nonisothermal  Batch  Reactor 

2.4.5  Construction  of  a  Digital  Controller 

2.4.6  Hybrid  Optimization  of  the  Sutro  Weir  Control  System 

2.4.7  Hybrid  Solution  of  Partial  Differential  Equations 

For  each  problem,  a  brief  introduction  is  presented  outlining 
the  nature  and  scope  of  the  simulation.  Next,  the  function  of  the 
analog  and  digital  portions  of  the  simulation  are  outlined.  Some 
examples  include  a  brief  description  of  the  theory  associated  with  the 
problem.  Finally,  simulation  results  are  presented  and  conclusions 
are  drawn.  More  information  on  the  simulation  programs  (flowcharts, 
analog  diagrams,  program  listings,  execution  instructions,  etc.)  may  be 
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found  in  a  research  report  which  describes  the  examples  in  detail  [38]. 

2.4.1  Digital  Initiation  and  Check  of  an  Analog  Simulation 
PURPOSE 

To  demonstrate  the  capability  of  a  hybrid  computer  and  speci¬ 
fically  the  Hytran  Operations  Interpreter  (HOI)  digital  language  in 
assisting  with  analog  computer  set-up  and  check  out.  The  analog  simula¬ 
tion  involves  the  two  tank  Sutro  Weir  experiment  used  in  the  under¬ 
graduate  laboratories. 

DESCRIPTION 

The  Sutro  Weir  experiment  is  shown  schematically  in  Figure  2.4 


FIGURE  2.4  SUTRO  WEIR  EXPERIMENT 


The  process  consists  of  a  single  feedback  control  system  used  to  control 
the  level  in  the  second  tank,  h^,  by  manipulating  the  flow  to  the  first 
tank,  q.  A  pneumatic  proportional  plus  integral  controller  is  used  in 
the  laboratory.  The  process  is  easily  modelled  with  simple  mass  balances 
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about  each  tank  and  controller  yielding  the  following  transfer  function 
relationships  between  levels  and  flow  [38]: 


H-j(s) 


M 


A-j  R-j  s+1 


+ 


hi 


A-j  R-j  s+1 


H2(s) 


H-|  (s)R-| 
A2R2S+I 


(2.4) 

(2-5) 


Q(s)  =  (R(s)  -  H2(s))(Kc  Kv(1  +  ^-)  (2.6) 

where  s  is  the  Laplace  transform  variable 

Q»L,H,R  are  variables  in  the  s-domain  as  described  by  Figure  2.4 
K  is  the  controller  gain,  K  is  valve  and  transmitter  gain 
Tj  is  integral  time*  A  is  the  area  of  a  tank 

Figure  2.5  illustrates  the  transfer  function  relationships  in 
block  diagram  form 


FIGURE  2.5  BLOCK  DIAGRAM  OF  TWO  TANK  CONTROL  SYSTEM 
ANALOG  FUNCTION 

To  simulate  the  Sutro  Weir  experiment  by  solving  the  dynamic 
equations  (2.4),  (2.5)  and  (2.6). 
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DIGITAL  FUNCTION 

To  use  the  HOI  language  to  calculate  and  set  potentiometers, 
to  make  a  static  check  of  the  analog  patch  panel,  and  to  make  a  dynamic 
check  of  the  analog  solution  of  the  differential  equations.  The  HOI 
language  is  an  interpretive  language  that  has  many  capabilities  for 
digitally  initiating  and  checking  analog  simulations  [5].  In  setting 
up  the  analog  simulation,  the  digital  program  calculated  and  set 
potentiometers  from  parameter  values  (controller  constants)  entered 
from  the  teletype.  This  saved  time  and  eliminated  human  error. 

The  static  check  was  automatically  performed  by  the  digital 
computer  (i.e.  HOI).  The  equations  describing  the  system  were  entered 
and  a  voltage  check  around  each  component  was  made  against  the  equations. 
This  effectively  checked  the  patching  and  the  components  in  the  analog 
simulation  [38]. 

The  dynamic  check  involved  a  numerical  solution  of  the  dif¬ 
ferential  equations.  The  equations  describing  the  Sutro  Weir  experiment 
used  in  the  static  check  were  solved  numerically  using  finite-difference 
approximations.  The  solution  was  converted  to  a  voltage  via  a  DAM  and 
plotted  on  the  analog  plotter  for  comparison  with  the  analog  results 
[38]. 

RESULTS 

The  results  of  the  analog  simulation  of  the  Sutro  Weir  are  pre¬ 
sented  in  Figure  2.6  for  load  and  setpoint  changes.  The  analog  solution 
was  obtained  after  a  static  check  of  the  analog  circuit  was  completed. 

The  static  check  involved  two  parts;  the  first  section  was  programmed  to 
describe  the  analog  patching  and  checked  the  analog  circuit  against  the 
dynamic  equations.  The  second  part  of  the  static  check  put  the  analog  in 
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FIGURE  2.6  SIMULATED  RESPONSE  OF  THE  SUTRO  WEIR  EXPERIMENT  TO  LOAD  AND  SETPOINT  CHANGES 
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the  Static  Check  mode  and  placed  voltages  at  the  input  to  components  in 
the  patch  panel;  and  checked  the  component  output  against  those  expected 
from  the  differential  equations  [38]. 

The  first  section  of  the  static  check  gave  no  errors  in  the 
analog  patching;  however,  the  second  section  revealed  two  errors.  A 
derivative  input  to  an  integrator  was  in  error  and  the  voltages  around 
a  potentiometer  failed  to  check.  The  potentiometer  setting  was  in 
error  and  was  manually  corrected.  The  balance  on  an  integrator  amplifier 
was  also  in  error  and  this  was  also  corrected. 

The  dynamic  check  of  the  response  is  presented  in  Figure  2.7. 
The  controller  constants  are  the  same  as  curve  1  in  Figure  2.6  for  a 
setpoint  change.  The  dynamic  check  uses  a  numerical  integration 
(trapazoidal  rule)  and  numerical  errors  could  arise.  The  calculation 
time  required  to  produce  this  solution  was  large  using  the  HOI  inter¬ 
preter  because  this  language  executes  source  coding.  As  a  result, 
execution  time  is  several  orders  of  magnitude  slower  than  compiler  or 
assembler  languages.  However,  the  source  code  is  easily  changed  on¬ 
line,  providing  very  simple  program  debugging.  Because  the  execution 
is  slow,  highly  accurate  numerical  solutions  are  not  feasible  for  the 
dynamic  check  of  large  problems.  For  large  problems  a  digital  machine 
using  an  analog  simulator  program  is  advisable.  To  illustrate  this,  a 
CSMP  (Continuous  Systems  Modelling  Program)  [39]  check  was  also  made 
and  the  result  is  shown  in  Figure  2.7.  The  CSMP  numerical  check  is  the 
more  accurate  for  it  uses  smaller  time  steps  and  a  more  accurate  second 
order  integration  approximation  than  the  HOI  check. 

CONCLUSIONS 


The  use  of  a  digital  computer  can  greatly  reduce  the  time 
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FIGURE  2.7 


HOI  AND  CSMP  DYNAMIC  CHECKS  OF  SUTRO  WEIR  SIMULATION 
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required  to  check  an  analog  simulation.  A  simple  interpretive  language 
is  desirable  since  it  is  easily  changed  and  high  speed  calculations  are 
not  required.  The  ability  of  the  language  to  provide  a  static  check 
provides  a  safeguard  against  patching  errors  and  component  failure.  A 
dynamic  check  gives  an  indication  of  the  correctness  of  the  final 
analog  solution. 

2.4.2  Reconstruction  of  Sampled  Signals 
PURPOSE 

To  demonstrate  the  errors  that  may  arise  in  sampling  and 
reconstructing  a  continuous  signal. 

DESCRIPTION 

A  sinusoidal  signal  is  generated  on  the  analog  computer.  It 
is  sampled  at  a  frequency  specified  by  the  user  and  stored  in  the  digital 
computer.  A  continuous  signal  is  then  reconstructed  from  the  stored 
values,  using  a  digital  to  analog  multiplier  (DAM)  and  its  associated 
zero  order  hold.  Finally,  a  sinusoid  of  the  same  frequency  is  calculated 
digitally  and  outputted  via  the  DAM.  Then  the  continuous  signal,  the 
sampled  signal,  and  the  digitally  generated  signal  can  be  compared 
graphically  from  the  X-Y  plotter  results. 

ANALOG  FUNCTION 

To  generate  a  sine  wave  of  amplitude  two  volts  and  frequency 
tt/4  radians  per  second  [38]. 

DIGITAL  FUNCTION 

To  acquire,  store,  and  playback  data  generated  on  the  analog 
at  various  sampling  frequencies  and  to  digitally  generate  A  sin  <*>t.  A 
FORTRAN  program  is  used  with  timing  controlled  by  core  cycle  time. 

The  digital  program  also  controls  the  analog  modes  necessary  to  plot  the 
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sinusoidal  functions  on  the  X-Y  plotter  with  the  DAM.  This  involves 
synchronizing  the  X-axis  sweep  (via  analog  integration)  with  the  start 
of  the  DAM  output  [38]. 

RESULTS 

Typical  results  are  illustrated  in  Figures  2.8  and  2.9.  As 
expected,  the  analog  sine  wave  is  more  accurately  reproduced  at  higher 
sampling  frequencies.  At  two  samples  per  second,  the  output  still 
gives  a  good  approximation  of  a  sine  wave.  Above  25  samples  per 
second  the  digitally  generated  function  appears  almost  continuous. 

Care  must  be  taken  to  ensure  that  the  sampling  frequency  is 
more  than  double  the  highest  frequency  of  the  sampled  signal  (Shannon's 
rule)  if  the  sampled  signal  is  to  be  recovered  completely.  For  example, 
this  sine  wave  would  have  to  be  sampled  at  least  once  every  four 
seconds  (natural  period  of  8  seconds)  to  recover  the  sine  wave  form. 

2.4.3  Simulation  of  Time  Delays 
PURPOSE 

To  demonstrate  the  simulation  of  time  delays  on  a  hybrid  com¬ 
puter. 

DESCRIPTION 

Time  delays  (or  transport  lags)  occur  frequently  in  process 
models.  Analog  simulation  of  transport  lags  requires  an  approximate 
representation  which  often  uses  many  operational  amplifiers  and  gives 
"good"  results  only  for  a  process  which  acts  as  a  low  pass  filter. 

Even  then,  undesirable  harmonics  caused  by  the  approximation  to  the 
time  delay  may  appear  in  the  simulated  process  response.  Hybrid  com¬ 
puters  can  use  the  digital  section  to  control  the  duration  of  time 
delays  precisely. 
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FIGURE  2.8 


SINE  WAVE,  25  SAMPLES  PER  SECOND 


33 


FIGURE  2.9  SINE  WAVE,  2  SAMPLES  PER  SECOND 
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Two  problems  of  this  type  are  demonstrated.  The  first 
problem  generates  a  time  delay  of  t  seconds  in  the  original  function, 

A  sinwt.  The  FORTRAN  subroutine  for  the  time  delay  was  written  to 
be  quite  general  and  can  be  applied  to  any  function,  f(t),  to  produce 
the  time-delayed  function,  f(t-x).  The  second  problem  simulates  the 
open  loop  dynamic  response  of  the  Sutro  Weir  experiment  for  several 
types  of  forcing  functions  using  the  HOI  language.  The  Sutro  Weir 
experiment  contains  a  fictitious  time  delay  causing  the  results  of  the 
forcing  function  to  be  seen  after  x  seconds.  Both  sections  are  com¬ 
pared  to  analog  simulated  time  delays. 

ANALOG  FUNCTION 

To  generate  a  sine  wave  and  to  simulate  the  Sutro  Weir  experi¬ 
ment  (without  controller)  for  open  loop  dynamic  testing  [38].  Timing 
for  the  transport  lag  in  the  Sutro  Weir  problem  is  controlled  by  inte¬ 
grator  monitoring  via  an  ADC  channel  [38],  and  has  a  resolution  of  .1 
second.  The  forcing  function  is  outputted  by  a  DAM. 

DIGITAL  FUNCTION 

To  provide  a  time  delay  in  a  sine  wave  and  the  Sutro  Weir 
simulation.  Timing  for  the  first  example  is  controlled  by  core  cycle 
time.  The  time  delay  subroutine  can  be  used  to  delay  any  analog 
function  to  within  t  .01  sec.  The  sine  wave  in  this  example  is  recon¬ 
structed  with  a  DAM  and  its  associated  zero  order  hold. 

RESULTS 

Figure  2.10  illustrates  the  hybrid  output  for  the  sine  wave 
of  example  2.4.2  with  a  time  delay  of  2  seconds.  Figure  2.11  illus¬ 
trates  the  same  function  using  CSDAP  [40]  (Control  System  Design  and 
Analysis  Program),  a  digital  continuous  process  simulator  that  uses  a 
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FIGURE  2.10  HYBRID  GENERATED  TIME  DELAYED  FUNCTION  Y  = 
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second  order  Pade  approximation  for  the  time  delay. 

Figure  2.12  illustrates  the  open  loop  response  of  the  Sutro 
Weir  with  a  two  minute  time  delay  to  step  and  pulse  changes  in  flow  to 
the  first  tank.  Figure  2.13  illustrates  the  same  response  using  CSDAP 
with  a  second  order  Pade  approximation  for  the  time  delay. 

CONCLUSIONS 

Both  hybrid  simulations  show  improved  time  delayed  responses 
over  the  analog  simulated  time  delays.  The  Pade  approximation  gives 
good  results  over  a  limited  frequency  range.  The  sine  wave  is  closer 
to  this  range  and  hence  the  adequate  results  shown  in  Figure  2.11. 

The  forcing  functions  of  the  second  program  have  a  greater  frequency 
content  than  the  sine  wave  and  exhibit  poorer  results.  The  high  fre¬ 
quency  ripples  caused  by  the  Pade  approximation  are  clearly  evident  in 
Figure  2.13. 

The  time  delay  simulated  in  the  HOI  language  is  a  special  type 
because  of  the  nature  of  the  forcing  functions  used.  These  functions 
do  not  change  with  time  after  they  are  inputted  to  the  system  and  can 
be  delayed  by  delaying  the  forcing  function  into  the  system.  In  the 
second  simulation,  the  time  delay  is  fictitious  and  has  no  physical 
equivalence  in  the  Sutro  Weir  equipment. 

The  FORTRAN  subroutine  used  in  the  first  example  is  more 
general  and  could  be  applied  to  plants  that  had  a  time  delay  with 
physical  significance.  For  example,  if  there  is  a  time  delay  between 
two  plants  whose  transfer  functions  are  known,  a  program  similar  to  the 
one  used  for  the  sine  wave  could  be  written  to  connect  the  two  plants, 
inputting  the  first  plant's  variable  to  the  second  plant  after  the 
desired  time  delay.  Reconstruction  of  time  delayed  signals  with  HOI  is 
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FIGURE  2.12  HYBRID  TIME  DELAYED  RESPONSE  OF  SUTRO  WEIR 
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FIGURE  2.13  CSDAP  TIME  DELAYED  RESPONSE  OF  SUTRO  WEIR, SECOND  ORDER 
PADE  APPROXIMATION 
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not  practical  due  to  slow  HOI  execution  time  and  present  HOI  timing 
considerations. 

2.4.4  Nonisothermal  Batch  Reactor 
PURPOSE 

To  compare  the  hybrid  and  CSMP  simulations  of  a  nonlinear 
batch  reactor. 

DESCRIPTION 

The  major  advantage  of  digital  simulation  over  analog  simula¬ 
tion  is  the  ability  of  the  digital  machine  to  handle  nonlinear  functions. 
The  errors  and  time  involved  in  numerical  integration  are  compensated 
for  by  the  accuracy  of  the  nonlinear  model  as  compared  to  a  linearized 
analog  model.  Nonlinear  analog  components  are  available  for  many 
functions,  however,  they  are  usually  less  accurate  and  more  expensive 
than  linear  components.  Hence  the  popularity  of  linear  analog  simula¬ 
tions.  A  hybrid  computer  system  has  the  ability  of  digitally  generate 
a  nonlinear  function  and  still  use  the  speed  of  the  analog  integrators. 

The  nonisothermal  batch  reactor  model  contains  an  exponential  nonlinearity 
since  the  assumed  reaction  rate  has  the  standard  Arrhenius -type  tempera¬ 
ture  dependence. 

ANALOG  FUNCTION 

To  simulate  a  nonlinear,  nonisothermal  batch  reactor  (see 
Derivation  of  Model ) . 

DIGITAL  FUNCTION 

To  provide  a  digital  simulation  of  the  nonlinear  batch  reaction 
using  CSMP  and  to  digitally  generate  the  function,  e^/RT),  for  use  in 
the  analog  simulation  [38].  An  analog  exponential  function  generator 
could  be  used  to  completely  solve  this  problem,  however,  many  functions 
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cannot  be  adequately  described  by  analog  components.  FORTRAN  is  used  as 
the  program  language  because  speed  is  essential  if  the  nonlinear  func¬ 
tion  which  is  calculated  and  outputted  to  the  analog,  is  to  closely 
approximate  the  correct  continuous  function.  HOI,  because  of  its  inter¬ 
pretive  nature,  has  a  slower  execution  speed  and  hence  the  results  using 
this  language  would  be  less  accurate. 

RESULTS 

Figure  2.14  compares  the  hybrid  simulation  for  reactor  tempera¬ 
ture  with  the  CSMP  simulation.  Figure  2.15  compares  hybrid  and  CSMP 
simulations  for  reactor  concentration.  The  computer  time  required  for 
10  hours  batch  reactor  simulation  time  is  30  to  282  seconds  (depending 
on  the  accuracy  desired)  using  CSMP  (IBM  1800  version)  and  10  seconds  on 
the  EAI  590  hybrid  computer. 

CONCLUSIONS 

Hybrid  simulation,  although  initially  more  time  consuming  to 
set  up  (patching,  magnitude  and  time  scaling),  can  have  substantial  time 
and  cost  savings  over  digital  simulation  if  multiple  runs  are  expected 
(e.g.  optimizing  batch  reactor  temperature  control).  Both  CSMP  and 
hybrid  simulations  can  easily  generate  arbitrary  nonlinear  functions. 
However,  the  time  involved  getting  accurate  numerical  results  (small 
integration  interval)  may  prohibit  many  digital  simulations. 

DERIVATION  OF  MODEL 

The  problem  chosen  for  simulation  is  an  ideal  first  order 
gas  reaction  carried  out  in  a  nonisothermal  batch  reactor. 
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TIME  IN  HOURS 


FIGURE  2.14  CSMP  VS  HYBRID  SIMULATION  OF  A  NONLINEAR  BATCH  REACTOR 
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FIGURE  2.15  CSMP  VS  HYBRID  SIMULATION  OF  A  NONLINEAR  BATCH  REACTOR 
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FIGURE  2.16  NONISOTHERMAL  BATCH  REACTOR 

Reaction  A  -*■  B  +  C 
Problem  Parameters: 

Rate  of  Reaction  rA  =  -  knA  =  -.856X10"^°  exp  (-1 .394X10^/T)hr~^n^  (2.7) 

Reactor  Volume  =  V  =  8.75  ft^ 

Cv>  Heat  Capacities:  A  =  30,  B  &  C  =  25  BTU./lb.mole  °F 
Heat  of  Reaction  =  AHrxn  =  -2500  BTU/(lb  mole  A) 

Initial  Temperature  =  600  °R 

Overall  Heat  Transfer  Coefficient  =  U  =  3.0  BTU/hr.ft^  °R 

Initial  Charge  =  nA  =  0.10  lb  moles 

2 

Heat  Transfer  Area  =  A  =  1  ft 
Material  Balance  on  A: 

dnA 

In  +  Generation  =  Out  +  Accumulation,  -  knA  =  (2.8) 

Energy  Balance: 

0  =  UA(T-Tm)  +  4Hrxn  (-kCAV)  +  Cv  £  (2.9) 

=  nA  CVA  +  nB  CVB  +  nC  CVC 

Substituting  the  parameters  gives  the  following  two  differential  equations: 
dnA  -in  4 

-^=  [-.856X10  lu  exp  (-1 .394X10  VT)]  nA  (2.10) 

JT  dn  n 

(5-20  nA)  ^  =  2500  -  3.0(T-540) 


(2.11) 
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2.4.5  Construction  of  a  Digital  Controller 
PURPOSE 
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To  program  a  general  digital  controller  on  a  hybrid  computer 
for  use  in  sampled-data  control  studies. 

DESCRIPTION 

A  digital  controller  is  constructed  for  the  analysis  of  control 
problems  which  result  when  the  controller  input  is  a  sampled  signal 
obtained  from  a  continuous  process.  The  analysis  of  such  sampled-data 
systems  is  usually  achieved  using  Z-transforms  [16].  In  this  example 
the  coefficients  of  the  Z-domain  digital  controller  are  specified  by  the 
user  and  are  inputted  to  the  digital  program.  The  controller  is  designed 
to  provide  a  "deadbeat"  system  response  [17]  for  a  first  order  plant, 
however,  the  controller  is  general  enough  for  use  in  other  types  of 
sampled-data  control  problems.  An  introduction  to  sampled-data  control 
systems  may  be  found  in  Chapter  13,  Automatic  Control  Engineering,  by 
F.  H.  Raven  [41]. 

ANALOG  FUNCTION 

The  continuous  first  order  plant  G(s)  =1/  (ts  +  1)  is  simulated 
on  the  analog  and  the  controlled  variable  is  sampled  by  an  ADC  channel. 
DIGITAL  FUNCTION 

To  act  as  a  digital  controller  (see  DIGITAL  CONTROLLER  THEORY 
section).  Direct  programming  [7]  was  used  in  the  controller  simulation. 
The  digital  computer  also  controls  the  sampling  of  the  plant  output, 
computes  the  control  action,  and  outputs  this  control  action  to  the 
plant. 

Both  HOI  and  FORTRAN  programs  are  used  with  FORTRAN  sampling 
being  controlled  digitally  by  core  cycle  time  with  resolution  to  1 
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millisecond.  The  HOI  sampling  is  controlled  by  integrator  monitoring 
with  resolution  to  0.1  seconds.  The  time  required  for  calculating  the 
control  action  upon  receiving  the  plant  error  and  implementing  the 
control  action  should  be  only  a  small  fraction  of  the  sampling  period. 
Again,  the  interpreter  (HOI)  execution  time  can  prove  prohibitive  for 
small  sampling  periods  [38]. 

RESULTS 


Results  for  the  first  order  plant  designed  for  a  deadbeat 
response  to  a  step  input  with  a  sampling  period  of  1  second  are  shown 
in  Figure  2.17.  This  figure  shows  the  different  responses  for  the 
FORTRAN  controller  and  the  HOI  controller.  Deadbeat  controllers  are 
extremely  sensitive  to  controller  parameters  and  timing  inaccuracies  [17], 
and  hence  the  slightly  inaccurate  deadbeat  response  found  in  the  HOI 
program  (due  to  timing  inaccuracies  and  significant  HOI  execution  time). 
The  HOI  program  could  be  used  in  less  critical  applications  with 
larger  sampling  periods. 

CONCLUSIONS 

Because  of  the  presence  of  the  zero-order  hold  in  the  DAM  cir¬ 
cuit  of  the  EAI  693  hybrid  interface,  the  hybrid  computer  is  ideally 
suited  for  the  study  of  sampled-data  control  problems.  Care  should  be 
taken  when  accurate  results  are  dependent  on  calculation  time  and  in 
these  instances  FORTRAN  Integer  Subroutines  or  Assembler  language  pro¬ 
grams  should  be  considered. 


DIGITAL  CONTROLLER  THEORY 


The  pulse  transfer  function  form  a  digital  controller  takes  is 
a„  +  a-,z”^  +  a0z*^  +  ...  +  a  z~n 

D(z)  =  - 1  2  " 


7]  I? 

+  bnz  +  b«z  +  . . .  +  bz 
o  1  2  m 


-m 


•  •  * 


(2.12) 
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where  b  4  0 
o 

The  control  system  is  represented  in  block  diagram  form  by  Figure  2.18 

/ 


I 

C(s) 


FIGURE  2.18  SAMPLED-DATA  CONTROL  SYSTEM 


D(z) 


E2(z)  _  aQ  +  a^"1  +  a2z'2  +  ... 

21 ^  b  +  b,z~^  +  b,z  2  +  . . . 
0  I  2 


(2.13) 


Cross  multiplying  yields 

b0  e2(z)  +  biz_1  E2(z)  +  b2z"2  E2^z^  +  +  bmz"m  E2(z)  = 

a0Ei  (z)  +  a^"1  E-j  (z)  +  a2z"2  E-j(z)  +  ...  +  anz"n  E^z)  (2.14) 
Taking  the  inverse  Z-transform  gives  the  time  domain  solution  as 
bQe2(t)  +  bi e2 ( t-T )  +  b2e2(t-2T)  +  ...  =  a  e^(t)  +  a^^t-T) 


+  . . .  +  a  e-j  (t-nT)  (2.15) 

or 

n  m 

4(t)  =  \ T  l  ake^t-kT)  ‘  IT  I  bke*(t-kT)  (2.16) 

2  bo  k=0  k  1  o  k=l  k  2 

■k 

where  e2(t)  is  the  control  action  at  the  present  sampling  instant 
e-|(t)  is  the  plant  error  at  the  present  sampling  instant 
T  is  the  sampling  period 
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2.4.6  Hybrid  Optimization  of  the  Sutro  Weir  Control  System 
PURPOSE 

An  optimization  problem  is  used  to  illustrate  a  class  of 
hybrid  computer  applications  in  which  multiple  analog  runs  are  performed 
and  the  digital  computer  is  used  both  to  monitor  the  analog  during  a  run 
and  perform  calculations  between  analog  runs.  Automatic  plotter  opera¬ 
tion  for  multiple  analog  runs  is  also  demonstrated. 

DESCRIPTION 

In  this  example,  "optimum"  controller  constants  are  calculated 
for  the  Sutro  Weir  control  system  described  in  2.4.1.  Optimization 
techniques  are  used  to  determine  controller  settings  for  a  P  +  I  con¬ 
troller  which  produce  "optimal"  process  responses  according  to  three 
different  performance  indices.  The  first  performance  index  is  minimum 
rise  time  subject  to  the  constraint  that  the  maximum  overshoot  be  less 
than  20%.  This  criterion  is  often  used  when  it  is  desired  to  drive  a 
process  to  within  specified  limits  in  the  shortest  possible  time.  A 
simple  exhaustive  search  is  carried  out  in  order  to  determine  the 
optimum  values  of  the  controller  constants,  Kc  and  ij. 

The  second  and  third  criteria,  the  integral  time  squared  error 
(ITSE)  and  the  integral  squared  error  (ISE),  are  common  in  many  control 
schemes.  A  systematic,  one  dimensional  search  is  used  to  discover  the 
values  of  the  proportional  gain,  K  ,  that  minimize  these  criteria. 

V 

ANALOG  FUNCTION 

To  simulate  the  Sutro  Weir  equipment  (Section  2.4.1)  with  the 
analog  output  h^  being  monitored  by  an  ADC  channel.  In  the  first  pro¬ 
gram  (minimum  time  criterion)  timing  is  controlled  by  monitoring  an 
integrator  via  an  ADC  channel.  In  the  second  program  (ISE  and  ITSE) 
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timing  is  controlled  by  a  "high"  signal  into  a  sense  line.  The  high 
signal  is  generated  by  the  analog  counter  [38]. 

DIGITAL  FUNCTION 

To  input  initial  controller  constants,  to  control  the  analog 
plotter  operation,  to  monitor  the  analog  during  a  simulation  run  in 
order  to  evaluate  a  performance  criterion,  to  calculate  new  controller 
constants  between  analog  runs,  and  to  use  the  calculated  controller 
constants  to  set  potentiometers  for  another  analog  run  [38]. 

RESULTS 

Figure  2.19  shows  the  output  of  one  of  the  six  runs  increment¬ 
ing  Tj  but  holding  KQ  constant  in  the  rise  time,  overshoot  optimization 
program.  From  the  results  of  the  runs,  a  value  of  K  =  2.75  and 

V 

tj  =  1.1  minutes  gives  a  minimum  rise  time  with  an  overshoot  of  18.75% 
[38].  Better  resolution  can  be  obtained  in  the  rise  time  calculations 
if  time  scaling  on  the  analog  circuit  is  decreased. 

Figure  2.20  shows  the  effect  of  Kc  on  the  ITSE  criterion  for  a 
step  change  in  setpoint  with  a  tj  value  of  1  minute.  A  final  Kc  value 
of  3.7  was  found  to  minimize  the  ITSE  criterion  for  a  setpoint  change 
of  6  inches.  Figure  2.21  shows  the  automatic  plotter  output  for  a 
load  change  of  0.7  cu.  ft.  per  minute.  A  controller  gain  of  14.5  (the 
gain  limit  of  the  physical  controller)  minimized  the  ITSE  criterion  for 
this  load  change  and  a  tj  of  1  minute. 

Figure  2.22  illustrates  the  plotter  output  for  the  minimization 
of  the  ISE  criterion  for  a  setpoint  change  of  6  inches  and  an  integral 
time  of  1  minute.  An  optimum  gain  of  6.9  was  discovered  in  this  case. 
CONCLUSIONS 

The  preceeding  results  illustrate  the  general  method  for  using 
the  hybrid  computer  to  solve  an  optimization  problem.  The  number  of  runs 
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TIME  IN  MINUTES 

FIGURE  2.19  OPTIMIZATION  OF  RISE  TIME  OVERSHOOT  CRITERION 
SETPOINT  CHANGE  .5  FEET 
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FIGURE  2.20  OPTIMIZATION  WITH  ITSE  CRITERION 
SETPOINT  CHANGE  OF  0.5  FEET  IN  H2 
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FIGURE  2.21  OPTIMIZATION  WITH  ITSE  CRITERION 

LOAD  CHANGE  OF  0.7  CU.  FEET  PER  MINUTE 
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FIGURE  2.22  OPTIMIZATION  WITH  ISE  CRITERION 
SETPOINT  CHANGE  OF  0.5  FEET 
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for  a  given  length  of  time  may  be  increased  by  operating  the  analog  in 
the  repetitive  operation  mode  or  by  increasing  integrator  speed.  Other, 
more  efficient,  optimization  techniques  (e.g.  golden  section  and  gradient) 
could  be  used;  however,  the  procedure  is  generally  the  same.  The  analog 
simulates  the  system  to  be  optimized,  while  the  digital  monitors  the 
analog,  evaluates  the  performance  criterion  and  makes  the  necessary 
optimization  calculations. 

2.4.7  Hybrid  Solution  of  Partial  Differential  Equations 
PURPOSE 

To  demonstrate  the  repetitive  use  of  an  analog  integrator  in 
the  hybrid  solution  of  partial  differential  equations  (P.D.E.). 

DESCRIPTION 

A  hybrid  computer  has  the  ability  to  use  an  analog  integrator 
repeatedly  in  order  that  problems  containing  many  integrations  may  be 
solved  with  a  limited  number  of  analog  components.  For  naturally  dis¬ 
crete,  time  varying  processes  (e.g.  plate  distillation  columns,  reactors 
in  series),  the  same  type  of  differential  equation  may  describe  the 
physical  situation  of  each  stage.  Therefore,  the  differential  equation 
describing  the  individual  stage  can  be  programmed  on  the  analog  and  the 
transfer  of  information  between  stages  can  be  controlled  by  the  digital. 
This  method  of  solving  repetitive  differential  equations  is  referred  to 
as  the  discrete-space-discrete-time  (DSDT)  technique  [6]. 

P.D.E.  may  be  solved  in  a  similar  manner  by  discretizing  one  of 
the  derivatives  (time  or  space)  and  solving  the  resulting  set  of  ordinary 
differential  equations  on  the  analog,  with  the  discrete  portion  solved 
digitally.  If  the  spatial  coordinate  is  solved  continuously  and  time  is 
incremented  digitally,  the  method  of  solution  is  referred  to  as  the 
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continuous-space-discrete-time  (CSDT)  technique.  The  CSDT  method  is 
illustrated  in  this  example  by  solving  the  P.D.E.  that  results  from  heat 
balances  on  a  tubular  heat  exchanger  (see  DERIVATION  OF  PROCESS  MODEL). 
ANALOG  FUNCTION 

To  solve  the  spatial  differential  equation  that  arises  when 
modelling  a  tubular  heat  exchanger  [38]. 

DIGITAL  FUNCTION 

To  increment  time  and  to  accumulate,  store  and  replay  data 
gathered  from  the  analog  in  order  to  numerically  solve  the  time  dis¬ 
cretized  differential  equation.  The  temperature  data  at  the 
sampling  instant,  T-|J  and  ,  are  sampled  every  0.01  feet  along  the 
length  of  the  exchanger  to  give  a  close  approximation  to  a  continuous 
function  in  the  spatial  direction  [38].  This  data  is  returned  to  the 
analog  as  and  l^"!  in  the  next  time  increment  (see  DERIVATION 

OF  PROCESS  MODEL). 

RESULTS 

The  temperature  profile  for  tube  fluid  temperature  is  illus¬ 
trated  in  Figure  2.23  (0.2  second  plotting  interval).  The  temperature 
profile  for  tube  wall  temperature  is  shown  in  Figure  2.24.  Figure  2.25 
illustrates  the  response  of  the  tube  fluid  temperature  to  a  step  change 
of  +100°F  in  the  tube  inlet  temperature. 

CONCLUSIONS 

Figure  2.23  compares  the  theoretical  steady  state  results  with 
the  hybrid  results  after  5  seconds.  The  agreement  is  excellent. 

The  solution  of  P.D.E.  on  a  hybrid  computer  has  been  studied 
in  great  detail.  At  very  rapid  solution  speeds,  sampling  error  can  be 
significant  and  many  published  papers  have  dealt  with  this  problem.  Also, 
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FIGURE  2.23  TUBE  FLUID  TEMPERATURE  PROFILE 


TIME  IN  SECONDS 
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FIGURE  2.24  TUBE  WALL  TEMPERATURE  PROFILE 
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TUBE  FLUID  TEMPERATURE  PROFILE, 
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difficulties  are  encountered  in  finding  a  value  of  the  time  increment 
At  that  yields  an  acceptable  compromise  between  the  effect  of  roundoff 
error  (in  analog  to  digital  and  digital  to  analog  conversion)  and  the 
truncation  error  associated  with  time  discretization.  However,  because 
of  the  overall  speed  and  accuracy  of  analog  integration,  hybrid  solu¬ 
tions  to  P.D.E.  are  usually  quite  accurate  and  can  be  obtained  much 
faster  than  by  pure  numerical  techniques. 

DERIVATION  OF  PROCESS  MODEL 

The  system  modelled  is  the  tubular  heat  exchanger  described  in 
Figure  2.26. 


v 

- > 

tube 
fluid  in 


v  7zzzzzzi£zz£zzzi2. 


I 

X 


X+AX 


r 


tube 

fluid  out 


FIGURE  2.26  SINGLE  PASS  HEAT  EXCHANGER 
An  energy  balance  on  the  increment  x  to  Ax  of  the  exchanger  yields 
(A,  v  +  h1P]  Axdg-Tp  =  (A,  v  PicpiTi)  Ix+Ax 

+  ft  <ixVlCplTl>  (2-17) 

Assume  the  physical  parameters;  A-j  =  area  of  tube,  p-j  =  density  of  tube 
fluid,  Cp-j  =  heat  capacity  of  tube  fluid,  P-j  =  perimeter  of  tube  available 
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for  heat  transfer;  are  constant.  Divide  by  ax  and  take  the  limit  as 
AX  -*  0 


Alplcpl  at1  +  Vlcplv  3X1  =  hlpl^T2  '  V 


(2.18) 


Rearranging 


3T,  3T, 

?r+vir  ■  ki^t2  -  Ti> 


(2.19) 


where 


K 


hipi 


i 


Vicpi 


A  similar  balance  on  the  tube  wall  assuming  T 2  to  be  independent  of  x 
gives 


aL 

3t~ 


=  K2CT3  -  T2)  -  K3(T2  -  T-j ) 


(2.20) 


where 


10 


h2P2 

A2p2Cp2 


and 


K 


hipi 


A2P  2Cp2 


T3  is  the  temperature  of  the  saturated  steam  and  therefore  is  constant 
with  time  and  length.  The  parameter  values  are  those  used  by  Vichnevetsky 
and  Tomalesky  [32]: 

L  =  10  ft  K-j  =  1  sec~^  K3  =  5  sec”^ 

v  =  2  ft/sec  K2  =  10  sec"^ 

A  numerical  backwards  difference  formula  is  used  to  approximate  the  time 
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derivatives.  Equations  (2.19)  and  (2.20)  become 


j  j  _  j  cjj  J 

At1  -  Mt23  -V)  (2-21) 

and 

j  j  _  j 

--At-2,1-  =  ¥T3°  '  T2°)  -  -  V>  <2-22> 

Rearranging  (2.21)  and  (2.22)gives  the  equations  solved  on  the  analog 


dT- 


dx 


K1  i  K1 

4  t,j  -  (4  + 


T  J’1 

_L)  T  j  +  h _ 

Atv'  1  Atv 


j  T2^”^  At  K2  T3j  K3At  . 

I  =  (!+AtK2+AtK3)  +  (l+AtK2+AtK3)  +  (l+AtK2+AtK3)  T1 


At  was  chosen  as  .05  seconds. 
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CHAPTER  THREE 

MODEL  REFERENCE  ADAPTIVE  CONTROL  BY  LIAPUNOV'S  DIRECT  METHOD: 

LITERATURE  SURVEY  AND  BASIC  EQUATIONS 

3.1  INTRODUCTION 

A  fixed  parameter  control  scheme  cannot  maintain  "good"  control 
in  industrial  processes  which  are  highly  nonlinear,  time  varying,  operate 
over  a  wide  range  of  process  specifications,  or  have  control  systems 
designed  on  the  basis  of  approximate  mathematical  models.  Consequently, 
frequent  tuning  of  control  loops  is  required  as  the  process  operating 
conditions  change.  Model  reference  adaptive  control  (MRAC)  is  a  control 
technique  that  changes  the  control  policy  to  partially  compensate  for 
varying  process  conditions.  A  reference  model  of  the  desired  closed-loop  pro¬ 
cess  dynamics  is  postulated;  the  adaptive  control  scheme  then  attempts 
to  minimize  the  error  between  the  process  response  and  the  reference 
model  response  when  both  are  subjected  to  the  same  inputs.  This  chapter 
presents  a  brief  historical  survey  of  MRAC  schemes,  from  single  variable 
gradient  techniques  to  the  multivariable  scheme  based  on  Liapunov 
stability  theorems  that  is  used  later  in  this  study.  The  theory  of 
multivariable  model  reference  adaptive  control  designed  by  Liapunov's 
direct  method  is  then  presented.  Extensions  to  existing  theory  are 
given  for  multivariable  adaptive  integral  and  setpoint  control.  Finally, 
the  theory  is  applied  to  the  design  of  a  MRAC  system  for  a  pilot  plant 
double  effect  evaporator. 

3.2  LITERATURE  SURVEY 

Model  reference  adaptive  control  systems  were  developed  in  the 
late  1 950 ' s  and  early  1 960 ' s  in  the  aircraft  and  aerospace  industries  for 
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applications  where  adequate  control  could  not  be  achieved  with  a  fixed 
parameter  control  design  [1].  In  many  applications,  because  of  model 
inaccuracy  and  process  nonlinearities,  the  control  scheme  had  to  change 
as  the  process  changed.  Eveleigh  ([1],  Chapter  13)  and  Bristol  [2] 
present  summaries  of  early  practical  applications  of  adaptive  control 
techniques  including  MRAC.  Bristol  also  outlines  the  early  single¬ 
input,  single-output  process  applications  using  conventional  analog 
controllers.  He  reported  that  adaptive  control  for  nonlinear  systems 
subjected  to  large  setpoint  changes  "works  in  a  trouble-free  and 
insensitive  manner  independent  of  process,  initial  settings,  distur¬ 
bances,  or  set-point  magnitude."  ([2],  page  561). 

The  early  model  reference  adaptive  control  systems  had  the 
same  objective  as  the  later  multivariable  techniques:  to  minimize  the 
error  (or  difference)  between  the  outputs  of  the  model  and  the  physical 
process  (see  Figure  3.1). 


FIGURE  3.1  BLOCK  DIAGRAM  OF  MODEL  REFERENCE  ADAPTIVE  CONTROL 
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The  early  techniques  used  sensitivity  functions  and  gradient-like 
searches  to  minimize  a  function  of  this  error  [1].  For  example  in  the 
"M.I.T.  technique"  [23],  the  design  objective  is  to  minimize  the 
integral  of  the  error  squared,  e  dt,  where  e  is  defined  as  the  error 

4 

between  the  model  response  and  the  process  response  for  the  same  input. 
A  differential  equation  for  adapting  the  controller  parameters  is 
obtained  based  on  the  gradient  of  the  error  with  respect  to  the 
adapting  parameters: 


da 

w  =  -  bey 

where  £  is  the  vector  of  parameters  being  adapted,  v  denotes  the 

_»Z_ 

gradient  with  respect  to  £,  and  b  is  a  positive  constant.  However, 
these  gradient  techniques  fail  to  guarantee  stability  for  even  simple 
systems  and  are  difficult  to  apply  to  multivariable  systems  [3]. 

Eveleigh  ([1],  Chapter  8)  describes  these  gradient  methods  in  greater 
detail  while  Parks  [3]  presents  informative  examples  showing  insta¬ 
bility  arising  from  adapting  with  the  M.I.T.  technique. 

In  an  attempt  to  overcome  the  stability  problems  associated 
with  gradient  techniques,  adaptive  schemes  based  on  Liapunov's  second 
or  direct  method  [4]  were  investigated  by  Grayson  [4],  Butchart  and 
Shackcloth  [5],  Parks  [3],  and  Shackcloth  [6]  for  simple,  single¬ 
variable  transfer  function  models.  Since  this  approach  is  based  on 
stability  considerations,  the  resulting  closed-loop  system  has  guaranteed 
stability  and  consequently,  higher  gains  can  be  used  in  the  adapting 
loops.  Simulation  studies  [5]  and  a  theoretical  analysis  [3]  demon¬ 
strated  that  the  Liapunov  approach  solves  the  stability  problems 
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encountered  using  the  gradient  methods.  Winsor  and  Roy  [7]  and  Porter 
and  Tatnall  [8]  extended  the  Liapunov  analysis  to  multivariable,  time- 
invariant  systems  presented  as  differential  equation  (state-space) 
models.  Porter  and  Tatnall 's  simulation  studies  [8]  demonstrated  a 
dramatic  improvement  over  fixed  parameter  control  schemes  with  the 
Liapunov  design  method  for  second  and  third  order  systems. 

The  "dual  problem"  of  process  identification  using  a  reference 
model  has  also  made  extensive  use  of  Liapunov  stability  considerations. 

In  the  process  identification  problem,  the  error  between  the  process 
and  the  reference  model  is  minimized  by  changing  the  model  parameters 
to  obtain  a  more  accurate  representation  of  the  physical  process. 
Apparently,  Rang  [9]  first  developed  the  Liapunov  model  reference  approach 
for  process  identification.  Grayson  [4]  presents  a  good  review  of  Rang's 
technique  and  related  methods.  Pazdera  and  Pottinger  [10]  demonstrated 
the  effectiveness  of  the  Liapunov  method  in  simulation  studies  and  were 
able  to  identify  the  process  parameters  using  only  process  noise  combined 
with  large  adaptive  loop  gains. 

More  recently,  Porter  and  Tatnall  [11]  have  extended  their 
adaptive  control  analysis  to  time-varying  systems  and  have  provided  an 
experimental  demonstration  of  this  technique  for  a  hydraulic  servo¬ 
mechanism  modelled  as  a  first-order  system  with  time-varying  gain  [12]. 
Their  experimental  results  tended  to  verify  the  simulation  studies  and 
showed  greatly  improved  control  over  a  fixed  parameter  scheme  because 
of  large  variations  in  process  parameters.  No  other  experimental  appli¬ 
cations  of  MRAC  systems  designed  by  Liapunov's  direct  method  have  been 
reported. 

Some  interesting  questions  related  to  MRAC  have  yet  to  be 
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answered.  One  of  these  is  discussed  by  Schooley  and  Pazdera  [13]  and 
Gromyko  and  Sankovskii  [24]  and  involves  the  proper  choice  of  the  adap¬ 
tive  loop  gains.  Since  the  controlled  system  is  guaranteed  to  be 
stable,  all  values  of  the  adaptive  loop  gains  can  be  used;  however,  in 
practice  some  stable  responses  are  more  desirable  than  others.  Schooley 
and  Pazdera  [13]  have  investigated  this  problem  and  have  attempted  to 
eliminate  loop  gain  sensitivity  by  adding  feedforward  control.  Gromyko 
and  Sankovskii  [24]  have  optimized  the  values  of  adaptive  loop  gain 
using  Pontryagin's  Maximum  Principle.  Their  work  illustrates  the 
improved  control  resulting  from  optimum  loop  gains  in  simulation  studies. 
This  problem  is  also  extensively  dealt  with  in  succeeding  chapters. 

A  few  applications  of  MRAC  to  process  control  problems  have 
been  reported  but  deal  exclusively  with  single-input,  single-output 
systems  which  are  presented  as  transfer  function  models,  controlled 
using  conventional  PID  controllers,  and  adapted  using  gradient  tech¬ 
niques.  Furthermore,  experimental  studies  have  only  been  reported  on 
the  earliest  and  simplest  techniques  [2].  Simulation  results  for  non¬ 
linear  continuous  stirred  tank  reactors  [14,15,16]  have  shown  great 
improvement  in  control  when  process  operating  conditions  vary  greatly 
(e.g.  start-up,  setpoint  changes);  however,  no  experimental  confirma¬ 
tion  of  the  simulation  results  has  been  reported. 

Crandall  and  Stevens  [14]  investigated  MRAC  of  a  single 
chemical  reaction  occurring  in  a  continuous  stirred  tank  reactor  (CSTR) . 
They  used  a  gradient  technique  to  adapt  the  gain  of  a  temperature  con¬ 
trol  loop;  their  simulation  studies  demonstrated  a  large  improvement 
over  a  fixed  gain  feedback  control  scheme.  Casciano  and  Staffin  [15] 
controlled  reactor  temperature  by  adapting  the  integral  time  of  a 


. 


■ 


68 


standard  proportional -integral  controller.  Again  improved  control  was 
reported.  Ahlgren  and  Stevens  [16]  applied  the  gradient  technique  to 
all  three  modes  of  a  PID  controller  for  the  same  reactor  example  used 
by  Crandall  and  Stevens.  They  also  simultaneously  adapted  the  para¬ 
meters  of  two  single  variable  controllers  for  reactor  temperature  and 
composition  respectively.  Excellent  control  was  achieved  over  a  wide 
range  of  process  operating  conditions.  They  also  demonstrated  their 
system  was  stable  for  certain  large  process  disturbances  but  did  not 
attempt  a  theoretical  investigation  of  system  stability. 

The  next  section  develops  a  MRAC  system  for  a  general  multi - 
variable  process  using  the  Liapunov  design  procedure  of  Porter  and 
Tatnall  [8]  to  guarantee  stability.  The  theory  is  then  modified  for  an 
actual  process  in  order  that  experimental  verification  of  this  method 
can  be  attempted. 

3.3  THEORY 

The  dynamic  behavior  of  many  plants,  including  the  double 
effect  evaporator  in  the  Department  of  Chemical  and  Petroleum  Engineer¬ 
ing,  can  be  adequately  approximated  by  a  state-space  model  of  the  form 

x  =  Ax  +  Bu  +  Dd  (3.1 ) 

"r  =  r  =  — 

where  Xp  represents  the  n  state  variables 
d  represents  the  p  disturbances 
u  represents  the  r  controls 

and  A,  B  and  D  are  constant  matrices  of  the  appropriate  dimensions. 

Assuming  multivariable  feedforward-feedback  control  of  the  form 


u 


Sfb  2p  +  £ff  - 


(3.2) 
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yields  a  closed-loop  process  model 

xp  =  (S  +  §  KFB)  *p  +  (g  +  |  KFF)  d  (3.3a) 


or  using  the  closed-loop  process  matrices  Ap  and  Dp  gives 


-P 


A_  xn  +  d. 

=p  -p  =p  - 


(3.3b) 


For  future  reference,  define  the  control  matrices  r  and  n  as 


l  -  SSfb 

n  =  B  Kpp 


(3.4) 


A  linear,  asymptotically  stable,  n^^-order  model  of  the 
desired  closed-loop  plant  behavior  is  postulated  as 


Am  x  +  Dm  d 
=m  -m  =m  - 


(3.5) 


where  xm  is  the  n-dimensional  state  vector  for  the  reference  model  and 
Am  and  Dm  are  specified  to  provide  the  desired  transient  and  steady 
state  behavior.  The  error,  e,  between  the  model  and  the  process  is 
defined  as 


The  intent  of  this  adaptive  control  technique  is  to  make  the  plant  and 
model  behavior  identical  by  changing  the  control  matrices  r  and  n.  The 
technique  for  changing  elements  in  r  and  u  develops  from  Liapunov's 
direct  method  of  stability  analysis  [4,8]. 

Subtracting  Equation  [3. 3b)  from  Equation  C3.5)  yields 
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S  =  te  +  ^  -  Sp)  Sp  +  (&■  -  2p)  4  <3-7) 

Consider  the  reduced  system: 

§  =  2m  e  (3.8) 

This  system  is  asymptotically  stable  since  the  model  in  Equation  (3.5) 
is  asymptotically  stable.  A  Liapunov  stability  theorem  states  that  for 
any  asymptotically  stable  linear  system,  such  as  that  described  by 
Equation  (3.8),  there  exists  a  positive-definite  matrix  P  that  is  the 
unique  solution  of 

aJ  p  +  p  A  =  -  g  (3.9) 

=m  =  =  =m  - 

where  Q  is  an  arbitrary  symmetric  positive-definite  matrix.  The 
Liapunov  function  is  then  chosen  as 

V  =  eT  P  e  (3.10) 

since  P  is  positive-definite.  The  total  time  derivative  of  the  Liapunov 
function  reduces  to 

V  =  -  eTQ  e  (3.11) 

which  is  negative-definite  since  Q  is  assumed  to  be  positive-definite. 
For  the  system  described  by  Equation  (3.7)  define 

Am  -  Ap  =  A  =  [c^j]  (i,j  =  l,2,...n)  (3.12) 

£m  '  5p  =  5  =  [X^j]  (1=1, 2, ...n;  j«l,2,...p)  (3.13) 


and 
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Following  Porter  and  Tatnall  [8],  choose  the  Liapunov  function  to  be  the 
positive-definite  quadratic  form 


n  n 


n 


v  =  eTPe  +  I  L'?L“,--i2+  I 


1*1  j=1  5ij 


1  x..2 


.  n  .L.  V  •  •  1  J 

1=1  J  =  1  1J  J 


(3.14) 


with  >  0,  v.  .  >  0;  £ .  .  and  v.  .  are  adaptive  loop  gains,  design 

'  J  *  J  »  J  l  J 

parameters  to  be  specified. 

The  Liapunov  function  defined  by  Equation  (3.14)  will  be  zero 
if,  and  only  if. 


e  =  0 


A  =  A 
=p  =m 


(3.15) 


D  =  D 
=p  =m 


The  total  time  derivative  of  Equation  (3.14)  is  given  by 


il  =  e^Pe  +  e’Pe  +  l  l 


n 


n 


n 

Z)  — —  a  •  •  a  •  •  +  )  )  -  A  .  .  A  .  . 

1  E •  •  lj  1J  •  n  V  •  •  1J  1J 

1=1  J=1  nj  J  J  1=1  J=1  1J  J  J 


A,,  A. 


(3.16) 


Substituting  Equation  (3.7)  for  e,  and  using  Equation  (3.9)  with 
Q  =  I  yields 


n  n 


v  = s  ?  +  V-i  j-Ei  +  xpJ-  0iJ 

n  p  T 

+  2  I  V  ( - A*.  +  d.e  P.)a.  • 

1=1  J=1  1J  J  J 


(3.17) 


th  * 

where  P^  denotes  the  i  column  of  matrix  P.  For  stability,  V  must  be 

negati ve-semi definite  which  requires  that  the  last  two  terms  of  Equation 

(3.17)  be  identically  equal  to  zero.  This  will  only  occur  if  a. .  =  0 

*  J 


' 
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and  a..  =  0  (i.e.  the  model  and  the  process  are  identical)  or  if 

*  vJ 


1  •  T 

a  •  •  +  X  .e  P  • 


1J  PJ- 


=  0 


(i  ,j  =  1  ,...,n)  (3.13) 


and 


1  •  T 

A.  .  +  d.e  P. 


vij  -1 


=  0  (i  =  1 , . .  ,n;  j  =  l,...,p)  (3.19) 


leaving 


V  =  - 


eT  e 


which  is  negative-semi  definite  in  e,  a.,  and  a...  Since  V  =  0  if,  and 

'  J  '  vJ 

only  if,  e  =  0  and  since  V  is  positive-definite,  it  follows  that  V  will 
decrease  monotonical ly  with  time  until  V  either  equals  zero  or  some 
positive  constant.  If  V  =  0  then  e  =  0  and  a.  .  =  0  and  A..  =  0;  alter- 

'  vJ  '  J 

natively,  if  V  eventually  attains  a  positive  value,  then  e  =  0  but 
a.| j  f  0  and/or  A^  f  0.  If  the  latter  case  occurs,  Am  -  Ap  =  Ca-jjl  1 
and/or  Dm  “  Dp  =  ^  0*  Consequently,  the  closed-loop  model  and 

process  matrices  do  not  have  to  be  equal  for  the  error  between  the  process 
and  model  states  to  be  zero.  This  point  is  emphasized  because  it  arises 
in  the  evaporator  application. 

Equations  (3.18)  and  (3.19)  can  be  rearranged  to  give 


a 


ij 


(3.20) 


and 


‘  dj 


Vj 


10 


(3.21) 


Equations  (3.20)  and  (3.21)  define  how  the  elements  of  the  matrices 
-m  ”  %  anc^  -m  "  -p  s^ou^  adapt.  However,  the  elements  which  can  be 
physically  manipulated  are  r  and  n  as  defined  by  Equation  (3.4),  The 
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corresponding  equations  for  the  elements  of  these  matrices  are  readily 
derived  by  differentiating  Equations  (3.12)  and  (3.13)  using  the  assump¬ 
tion  that  the  open-loop  plant  and  model  matrices  are  time-invariant. 


• 

Caij]  = 

^m  ’  =p 

-ft.  - 

-  (A  +  r)  = 

i — i 

•  >- 
i _ i 

• 

H 

•  p-  M 

i 

(3.22) 

^  = 

•  • 

D-D  = 
=m  =p 

• 

-=DP  ■ 

•  • 

-  (D  +  n)  = 

•  • 

-  U  =  -  [to.j  j  ] 

(3.23) 

Therefore,  Equations  (3.22)  and  (3.23)  combined  with  Equations  (3.20) 
and  (3.21)  yield  the  following  two  adaptive  control  laws  for  changing 
the  elements  of  r  and  n. 


(3.24) 


d.  e"*"p. 
J  -  Li 


V. 


1J 


(3.25) 


The  block  diagram  in  Figure  3.2  reviews  the  interrelations  of  the  equa¬ 
tions  in  this  adaptive  control  scheme. 


FIGURE  3.2  BLOCK  DIAGRAM  OF  MRAC  VIA  LIAPUNOV'S  DIRECT  METHOD 
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The  MRAC  scheme  based  on  Liapunov's  direct  method  can  be  con¬ 
sidered  to  be  nonlinear  due  to  the  adapting  algorithms  given  by  Equations 
(3.24) and  (3.25).  The  method  requires  measurements  or  estimates  of  all 
the  process  states,  Xp,  and  disturbances,  and  also  the  ability  to  change 
the  elements  of  r(t)  and  n(t). 

3.4  EXTENSIONS  TO  THE  EXISTING  THEORY:  ADAPTIVE  INTEGRAL  AND  SETPOINT 

CONTROL 

The  control  theory  just  presented  considers  feedforward  and 
proportional  feedback  control.  For  many  processes,  integral  control 
action  is  desirable  in  order  that  offset  in  some  of  the  variables  may  be 
eliminated.  Also,  setpoint  changes  cause  the  process  conditions  to 
change  and  thus  require  the  parameters  in  the  linear  process  model  to 
be  updated.  The  theory  for  adding  adaptive  integral  and  setpoint  con¬ 
trol  is  now  developed  as  an  extension  of  the  feedforward-feedback  theory 
just  presented. 

3.4.1  Adaptive  Integral  Control 

The  open-loop  state-space  model  of  the  process  is  again  des¬ 
cribed  by  Equation  (3.1) 

x  (t)  =  A  x  (t)  +  B  u(t)  +  D  d(t)  (3.1) 

It  is  desirable  to  eliminate  the  offset  in  some  of  the  state  variables, 
y,  a  subset  of  x  .  Let  the  dimension  of  y  be  q.  Newell  [17]  has  demon- 

r 

strated  that  for  all  q  "integral  states"  to  be  equal  to  zero,  q  must  be 
less  than  or  equal  to  r,  the  dimension  of  the  control  vector,  u.  We 
can  write 

t  =  £  *p 


(3.26) 


■  - 
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where  C  is  specified  to  make  y  the  desired  subset  of  xp. 
Define  a  new  set  of  state  variables  z  as 


z 


•t 

y  dt 
o 


or 


(3.27) 


1 


=  2 


(3.28) 


Equations  (3.28)  and  (3.26)  give 


l  =  £  xp 


(3.29) 


Augment  the  state  vector,  x^  with  the  new  states  z  by  using  Equations 

(3.29)  and  (3.1) 


• 

*p 

z 

= 

A  0 

C  0 

Sp 

z 

+ 

B 

0 

u  + 

D 

0 

zz  zz 

=: 

Define 


-P 

B 

D 

z 

;  A'  = 

C  0 

;  B'  = 

0 

ii 

Q  II 

•  r* 

1 

1  O  II  | 

(3.30) 


(3.31) 


A  feedback  control  law  is  now  used  for  the  augmented  system 


u(t)  =  K'  xj,(t)  (3.32) 

Define 

K'  =  [KpB,  Kr]  (3.33) 

Expanding  Equation  (3.32)  using  Equations  (3.31),  (3.28)  and  (3.27) 
yields  the  proportional  integral  control  action 
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rt 

u  =  KrR  x  +  Kt  y  dt  (3.34) 

=  FB  “P  =I  J0 

Substituting  Equation  (3.32)  into  Equation  (3.30)  yields 

Xp  =  [A‘  +  B'  Kl]  x’  +  D1  d  (3.35) 

Equation  (3.35)  is  now  in  the  form  of  Equation  (3.3a)  (without  feed¬ 
forward  control)  and  the  adaptive  control  law  defined  by  Equation  (3.24) 
can  be  used  to  generate  the  elements  of  K'  and  hence  Kp^  and  Kj . 

The  addition  of  integral  control  adds  more  elements  to  be  cal¬ 
culated  by  the  adaptive  control  law.  The  number  of  adaptive  control 
calculations  increases  as  the  square  of  the  augmented  state  vector  so 
integral  action  can  greatly  increase  the  adaptive  control  calculation 
time.  The  dimension  of  K1  is  n+q  X  r  instead  of  n  X  r  for  Kpg.  The 
matrices  for  the  augmented  reference  model 


(3.36) 


can  be  obtained  from  optimal  control  applied  to  the  augmented  state 
vector  [17]  or  by  other  techniques. 


3.4.2  Adaptive  Setpoint  Control 

Setpoint  changes  are  an  important  adaptive  control  applica¬ 
tion  because  when  a  setpoint  change  is  made  to  a  nonlinear  system,  the 
control  action  should  adapt  in  order  to  maintain  the  best  control  at 
the  new  process  operating  conditions.  A  control  law  that  has  been 
derived  from  an  open-loop  linear  model  at  one  set  of  operating  condi¬ 
tions  cannot  be  expected  to  work,  as  well  at  conditions  far  from  where 
the  linear  model  was  obtained. 
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The  open-loop  dynamics  of  the  system  can  again  be  described  by 
Equation  (3.1) 


(3.1) 


xn=  Ax+Bu+Dd 
-P  «  -p  =  -  =  - 


Let  the  state  variables  that  are  subjected  to  setpoint  changes  be 
defined  by 


(3.37) 


=  -P 


y 


Again  y  is  a  subset  of  >^as  defined  by  matrix  C.  The  maximum  dimension 
of  y  is  r,  the  dimension  of  the  controls  [17,25].  In  other  words,  only 
a  subset  of  the  state  variables  equal  to  the  number  of  control  variables 
can  have  their  setpoints  arbitrarily  specified  and  have  the  remaining 
states  stay  at  their  original  setpoint. 


A  control  law  that  accounts  for  setpoint  changes  is  of  the 


form 


y(t)  =  =fb  $,(t)  +  £ff  $(*)  +  =KsP  ysp(t) 


(3.38) 


where  y  is  the  desired  value  of  y  and  can  be  time-varying.  Using  the 
-sp  - 

setpoint  vector  to  augment  the  disturbance  vector 


u(t)  =  Kfb  x(t)  +  Kpp'  d ' (t ) 


(3.39) 


where 


(3.40) 


d 


(3.41) 


d 


Substituting  Equation  (3.38)  into  Equation  (3.1)  yields 


. 


. 


. 
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xp  =  (A  +  B  KFB)Xp+  (D  +  B  Kpp)d  +  B  £sp  ysp  (3.42) 

Defi  ne 

D'  =  [D,0]  (3.43) 

Substituting  Equations  (3.43),  (3.41)  and  (3.40)  into  Equation  (3.42) 
yields 


xp  =  (A  +  B  Kfb)^+  (D‘  +  B  Kpp 1 )  d'  (3.44) 

This  is  in  the  standard  form  of  Equation  (3.3a)  and  the  adaptive 
control  law  for  adapting  the  feedforward  matrix  as  described  by  Equation 

(3.25)  can  be  used.  The  resulting  control  configuration,  feedback 
feedforward  setpoint  control,  again  has  the  disadvantage  of  requiring 
additional  adaptive  control  calculations.  The  increase  in  time  required 
is  proportional  to  the  square  of  the  dimension  of  the  y  vector. 

The  model  matrices  in  the  closed-loop  model  equation 

x  =  Am  x+  D '  d '  (3.45) 

can  be  determined  from  a  modification  of  the  optimal  control  problem 

[17.25]  or  from  other  considerations  such  as  state  driving. 

Both  integral  and  setpoint  adaptive  control  can  be  used 
simultaneously  by  augmenting  the  state  vector  for  integral  control  and 
the  disturbance  vector  for  setpoint  control.  Newell  [17,25]  has  con¬ 
sidered  these  extensions  to  the  standard  optimal  regulator  problem. 

The  resulting  control  scheme  contains  feedforward  and  setpoint  modes 
as  well  as  proportional  plus  integral  feedback  control.  That  is  a 
control  law  results  of  the  form 


. 


* 
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u(t)  =  Kfb  j^t)  +  Kj 


y(t)dt  +  Kff  d(t)  +  Ksp  y$p(t) 


(3.46) 


3.5  APPLICATION  OF  MRAC  TO  A  PILOT  PLANT  EVAPORATOR 

This  section  demonstrates  how  the  theory  described  in  the 
previous  section  can  be  applied  to  the  double  effect  evaporator  in  the 
Department  of  Chemical  and  Petroleum  Engineering  of  the  University  of 
Alberta.  A  schematic  diagram  of  the  pilot  plant  evaporator  is  shown  in 
Figure  3.3.  The  first  effect  is  a  calandria  type  unit  with  32  eighteen 
inch  long,  3/4  inch  O.D.  tubes.  The  feed  to  this  three  gallon  unit  is 
5  pounds  per  minute  of  3  percent  by  weight  triethylene  glycol  and  this 
feed  is  heated  by  2  pounds  per  minute  of  250°F  saturated  steam.  The 
second  effect  has  externally  forces  circulation  through  three  6  feet 
long,  one  inch  O.D.  tubes  fed  by  the  first  effect  product  and  heated  by 
the  first  effect  vapour.  The  second  effect  vapour  is  totally  condensed. 
Vacuum  control  maintains  the  necessary  pressure  differential  between 
effects.  The  product  is  about  1.5  pounds  per  minute  of  10%  glycol.  A 
more  detailed  description  of  the  evaporator  can  be  found  in  a  thesis  by 
Newell  [17,  Chapter  4]  or  a  paper  by  Andre  and  Ritter  [19]. 

The  evaporator  has  been  extensively  modelled  by  Newell  [17], 
Wilson  [18],  and  Andre  and  Ritter  [19].  The  model  used  here  is  due  to 
Wilson  [18]  and  is  a  five  state,  linear,  time-invariant  model  in  the 
form  of  Equation  (3.1)  or 

xn  =  Ax+Bu+Dd 
-p  =  -p 

The  following  definitions  for  the  state,  disturbance,  and  control  vectors 
x  ,  d,  and  u,  were  used: 
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FIGURE  3.3  SCHEMATIC  DIAGRAM  OF  DOUBLE  EFFECT  EVAPORATOR 
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where  W-j ,  C-j ,  S',  F'  etc.  denote  normalized  perturbated  variables 


(wrwiss) 

i.e.  W-j  =  — Tj -  (ss  -  steady  state) 

W1  ss 


and  the  symbols  are  defined  in  the  Nomenclature  section  associated  with 
this  Chapter. 


Normal  process  operating  conditions  [18]  (presented  in  the 
Appendix  for  Chapter  Four)  gave  the  following  numerical  values  for  the 


A,  B  and  C  matrices. 


A 


0 

0 

0 

0 

0 


-.0011  -.1255 

-.0755  .1255 

-.0060  -.7741 

-.0012  -.1448 

.0393  .1448 


0  0 

0  0 

0  0 

0  .0001 

0  -.0380 


. 

■ 


' 


82 


0  -.0766  0 


B 


0  0  0 

.2160  0  0 


0  .0795  -.0381 

0  -.0414  0 


D 


.1098  0 

-.0333  .0766 

-.0188  0 

0  0 

0  0 


.0911 


0 

0 


0 

0 


For  implementation  on  the  evaporator,  only  multivariable 


proportional  feedback  control  was  considered  since  previous  studies  [17] 
have  demonstrated  that  this  configuration  gives  excellent  control. 
Therefore,  in  Equations  (3.3)  and  (3.4),  Kpp  =  0  and  n  =  0.  The  r  matrix 
generated  by  the  adaptive  control  law  in  Equation  (3.24)  consists  of  25 
elements.  However,  only  15  elements  of  KpB  can  be  adapted  since  B  is 
fixed  by  the  physical  process  (note  that  r  =  B  KpB  in  Equation  (3.4)). 
Consequently,  the  pseudo-inverse  approach  [20]  was  used  to  generate  the 
following  "solution"  to  Equation  (3.4); 


TD\-1  DT 


KFB(t)  =  (B'B)"'  B'  r(t) 


(3.47) 


Equation  (3.47)  provides  a  least  squares  solution  for  kpg(t) 
that  is  only  approximate  since  the  15  unknown  elements  of  Kpg  are 
generated  from  25  scalar  equations  involving  the  elements  of  r(t)  and  B. 

The  model  of  the  desired  closed-loop  behavior  is  chosen  to  be 
the  open-loop  model  plus  the  optimal  feedback  control  law  for  a  quad¬ 
ratic  performance  index.  This  optimal  feedback  matrix  is  derived  using 


. 
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dynamic  programming  and  minimizes  the  quadratic  performance  index  used 
by  Newell  [17,  Chapter  5]. 

The  r(t)  matrix  is  calculated  from  Equation  (3.24)  using  a 
backwards  difference  approximation.  This  yields 

Yij*(k)  =  xpj.(k)  ?(k)T  ^ij  At  +  YijCk-D  (3.48) 

The  initial  conditions  for  Equation  (3.48),  y.  .(0),  are  related  to  the 
initial  control  matrix,  Kpg(0),  and  must  be  specified.  The  importance 
of  the  initial  control  is  investigated  in  Chapters  Four  and  Five. 

The  values  of  x  are  sampled  from  the  continuous  process  and 

r 

the  model  states  are  calculated  from 

xm(k+1)  =  xm(k)  +  em  d(k)  (3.49) 

-m  im  -nr  1  =m  -v 

which  is  Equation  (3.5)  in  discrete  form.  The  numerical  values  of  the 


discrete  matrices 

^m 

and  are: 

=m 

"  .4341 

i 

• 

o 

1 

o 

-.0785 

0 

-.7874 

.1248 

.9037 

.0313 

0 

-.2266 

tm 

= 

1.692 

-.2564 

-.3186 

0 

-3.071 

0 

0 

0 

-.0001 

-.0003 

-.0835 

.0015 

.0150 

0 

.1515 

“  .1182 

0 

-.0050" 

-.0351 

.0785 

.0049 

E 

CD  II 

= 

-.0136 

-.0002 

.0662 

.0012 

0 

-.0058 

-.0019 

.0016 

.0058 

Details  on  the  derivation  of  these  matrices  are  given  in  the  Appendix  for 
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this  chapter. 

A  review  of  the  theory  presented  above  suggests  the  following 
procedure  for  the  implementation  of  multivariable  model  reference 
adaptive  control.  At  each  sampling  instant: 

a)  Sample  the  plant  state,  xp(k),  from  the  continuous  process. 

b)  Calculate  the  model  state  for  the  optimal  feedback  model, 
xm(k),  from  the  discrete  closed-loop  model  in  Equation 
(3.49). 

c)  Compare  these  states  and  from  the  error,  e(k)  =  xm(k)  - 
Xp(k),  calculate  the  control  matrices  r(k)  and  KpB(k) 
from  Equations  (3.47)  and  (3.48). 

d)  Calculate  the  control  variables  u(k)  from  u(k)  =  KpB(k) 
Xp(k)  and  send  these  signals  to  the  process. 
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CHAPTER  FOUR 

HYBRID  SIMULATION  OF  A  MODEL  REFERENCE  ADAPTIVE 
CONTROL  SYSTEM  FOR  THE  EVAPORATOR 

4.1  INTRODUCTION 

The  model  reference  adaptive  control  (MRAC)  scheme  described 
in  the  previous  chapter  was  simulated  on  the  EAI  590  hybrid  computer 
system  in  the  Department  of  Chemical  and  Petroleum  Engineering  at  the 
University  of  Alberta.  The  EAI  580  analog  computer  was  used  to  simulate 
the  dynamic  behavior  of  the  double  effect  evaporator,  as  represented  by 
the  fifth-order  linear  evaporator  model  of  Wilson  [2].  The  reference 
model  response  and  the  adaptive  control  matrices  were  generated  on  the 
EAI  640  digital  computer.  The  hybrid  simulation  was  ideal  for  evaluat¬ 
ing  this  computer  control  technique  because  the  analog  model  provided 
an  intuitive  feel  for  the  behavior  of  the  actual  evaporator,  while  the 
digital  computer  provided  an  authentic  simulation  of  the  process  control 
computer  required  for  actual  implementation.  Furthermore,  the  digital 
control  programs  used  in  the  simulation  were  easily  transferred  to  the 
IBM  1800  computer  for  experimental  verification  of  the  technique. 

4.2  ANALOG  EVAPORATOR  MODEL 

The  analog  simulation  of  the  evaporator  was  essentially  the 
same  as  the  one  used  by  Newell  [7]  with  modifications  in  the  disturbance, 
control,  and  data  monitoring  sections  for  control  by  the  EAI  640  digital 
computer.  The  state-space  evaporator  model  chosen  was  in  normalized 
perturbation  form  making  all  the  states,  controls,  and  disturbances  vary 
about  the  steady  state  value  of  zero  and  generally  stay  within  the  range 
of  t  1.0  [3].  This  model  was  ideal  for  analog  simulation  because 
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magnitude  scaling  was  not  required  unless  extremely  large  disturbances 
or  inadequate  control  schemes  were  to  be  considered.  Figure  4.1  shows 
the  analog  circuit  used  for  the  evaporator  simulation.  The  nonlinear 
evaporator  model  [5]  would  have  provided  more  accurate  process  dynamics; 
however,  a  lack  of  nonlinear  elements  prevented  this  simulation. 

Communication  with  the  digital  computer  was  achieved  through 
analog  to  digital  converters  (ADCs)  to  access  the  process  states  as 
required  by  the  MRAC  algorithm.  The  analog  simulation  was  time  scaled 
by  a  factor  of  60,  i.e.  one  second  of  simulation  time  represented  one 
minute  of  real  time.  Also,  data  were  acquired  automatically  for  later 
analysis  and  plotting  by  the  IBM  1800  computer  for  5  minutes  of  simula¬ 
tion  time.  More  details  on  the  analog  simulation  and  the  data  accumula¬ 
tion  programs  can  be  found  in  a  hybrid  user's  manual  by  Oliver  and 
Seborg  [4]. 

4.3  DIGITAL  COMPUTER  CONTROL  PROGRAMS 

The  digital  computer  program  controlled  sampling  of  the  evap¬ 
orator  states  from  the  analog  and  calculated  the  model  states  from  the 
discrete  form  of  the  reference  model  (Chapter  Three).  These  states  were 
compared  and  the  resulting  errors  and  acquired  process  states  were  used 
to  generate  the  adaptive  control  matrix,  r(t).  This  matrix  was  then 
used  to  calculate  the  feedback  matrix,  KpB,  and  subsequently  the  controls, 
u(t).  These  controls  were  then  sent  to  the  analog  evaporator  model  via 
digital  to  analog  multipliers  (DAMs).  The  logic  flowchart  for  these 
calculations  and  complete  program  listings  can  be  found  in  a  research 
report  [4]. 


4.4  RESULTS 
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FIGURE  4.1  ANALOG  SIMULATION  OF  THE  LINEAR,  FIFTH-ORDER  EVAPORATOR  MODEL 
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Simulation  runs  were  made  to  investigate  several  factors  which 
affect  the  performance  of  the  adaptive  control  system  including: 
initial  control  policy,  adaptive  loop  gains,  number  of  elements  of  r 
which  adapt,  choice  of  reference  model,  and  "unseen"  disturbances. 

Unless  otherwise  noted,  all  simulation  runs  considered  a  step  disturbance 
of  +  20%  in  feed  flow  and  assumed  that  initially,  the  evaporator  was 
operated  in  the  open-loop  configuration  (i.e.  no  control,  r(0)  =  0). 

The  horizontal  arrows  on  the  figures  indicate  the  desired  setpoint  of 
the  system;  the  vertical  arrows  represent  the  start  of  the  disturbance. 
The  most  noteworthy  results  were  as  follows: 

4.4.1  Improvement  in  Control 

Figure  4.2  illustrates  the  improved  control  resulting  from 
MRAC  and  displays  the  reference  model  response,  the  adaptive  control 
response,  and  the  open-loop  response  with  no  control  of  and  only 
proportional  control  of  the  two  levels.  The  adaptive  response  starts 
with  no  control  initially  (i.e.  r(0)  =  0),  a  policy  which  is  even  worse 
than  the  level  control  presented  in  Figure  4.2.  If  this  initial  control 
policy  were  continued  (without  adaptive  control),  the  first  effect  level, 
W-j ,  would  drop  to  zero  in  approximately  45  minutes.  Despite  this 
unfavorable  starting  point,  the  adaptive  control  scheme  quickly  brought 
the  process  response  back  to  the  model  response.  The  model  response  is 
excellent  since  it  utilizes  optimal  feedback  control  based  on  a  quadratic 
performance  index  [5]. 

4.4.2  Effect  of  on  MRAC 

—Li 

The  effect  of  the  adaptive  loop  gains,  is  demonstrated  in 

* 

Figure  4.3  where  values  of  1 ,  10  and  100  for  all  elements  of  5  were 
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FIGURE  4.3  SIMULATED  RESULTS  DEMONSTRATING  THE  EFFECT  OF  ADAPTIVE 
LOOP  GAIN,  ON  MRAC 
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investigated.  All  three  responses  are  stable  as  theory  predicts; 
however,  when  =  1,  control  appears  sluggish  and  for  £..  =  100,  the 

*  J *  *  J 

system  responds  in  an  oscillatory  manner.  When  £• .  .  =  10,  the  process 

’  J 

response  adapts  to  the  model  response  most  rapidly.  Thus,  as  in  con¬ 
ventional  feedback  control  systems,  there  appears  to  be  a  range  of 
gains  which  provide  a  satisfactory  trade-off  between  speed  of  response 
and  degree  of  oscillation.  On  the  basis  of  the  results  shown  in 
Figure  4.3,  a  gain  somewhere  between  10  and  100  would  be  expected  to 
give  the  best  response.  However,  no  attempt  was  made  to  find  the 
optimum  value  of  £  • . .  Figures  4.4  and  4.5  demonstrate  the  effect  of 

*  vJ 

£.jj  on  e-j  and  y^  and  illustrate  the  rapid  reduction  in  error  for 
=  10  and  the  more  sluggish  adapting  that  occurs  when  =  1. 

1  J  '  J 

To  demonstrate  that  the  controlled  system  remains  stable 
despite  large  values  of  adaptive  loop  gains,  5. .  was  set  equal  to  10,000 

■  J 

in  the  run  shown  in  Figure  4.6.  As  this  figure  demonstrates,  the  control 
is  quite  oscillatory  but  still  stable. 

4.4.3  Number  of  Adapting  Elements  of  r 

Figure  4.7  illustrates  that  all  elements  of  r  need  not  adapt 
for  MRAC  to  perform  satisfactorily.  The  curves  in  Figure  4.7  are  for 
successively  fewer  elements  of  r  adapting  (21,  11  and  7)  and  still  satis¬ 
factory  control  is  maintained.  However,  with  fewer  elements  adapting 
(i.e.  some  £ . .  =  0),  successively  worse  model  following  results.  Also, 

*  J 

it  is  increasingly  difficult  to  predict  what  will  happen  when  a  disturb¬ 
ance  occurs  (note  that  the  initial  response  changes  direction  when 
the  number  of  elements  adapting  is  changed  from  11  to  7). 

Figure  4.8  demonstrates  that  as  few  as  4  elements  need  adapt 
for  adequate  control  to  be  maintained.  However,  the  elements  of  r  which 
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FIGURE  4.4  SIMULATED  EFFECT  OF  5...  ON  e1  =  xml  -  xpl 
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FIGURE  4.5  SIMULATED  EFFECT  OF  5tJ  ON  RATE  OF  ADAPTING  FOR  y5] 
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FIGURE  4.6  SIMULATED  RESULTS  DEMONSTRATING  INHERENT  STABILITY 
OF  MRAC  DESIGNED  VIA  LIAPUNOV'S  DIRECT  METHOD 
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FIGURE  4.7  SIMULATED  RESULTS  FOR  ONLY  SOME  ELEMENTS  OF  r 
ADAPTING,  5..  =10 
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FIGURE  4.8  SIMULATED  RESULTS,  FOUR  ELEMENTS  OF  r  ADAPTING, 
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do  adapt  cannot  be  arbitrarily  chosen  because  with  some  £ .  .  =  0, 

*  vJ 

stability  is  no  longer  guaranteed.  For  example,  it  was  observed  that 
with  6  adapting  r  elements,  an  unfortunate  choice  of  adapting  elements 
would  cause  an  unstable  response.  This  situation  occurs  because  the 
open-loop  system  is  unstable  and  at  least  one  element  in  both  the  first 
and  fourth  rows  of  r  must  adapt  if  there  is  to  be  level  control  for  the 
first  and  second  effects.  The  successive  elimination  of  adapting  r 
elements  was  chosen  so  as  to  retain  the  subsequent  adapting  of  the 
largest  element  of  the  Kpg  matrix,  as  determined  from  examining  the 
optimal  feedback  matrix  used  in  the  reference  model.  The  Appendix 
associated  with  this  chapter  presents  a  table  showing  the  adapting 
elements  for  each  run. 

Figures  4.8  and  4.9  illustrate  that  higher  adaptive  loop  gains 
are  possible  without  causing  oscillatory  responses  when  some  adapting 
elements  are  eliminated.  Figure  4.9  shows  improvement  in  control  when 
£..  changes  from  10  to  100  for  7  elements  adapting.  By  contrast,  a  gain 

*  J 

of  100  produced  an  oscillatory  response  when  all  the  elements  of  r 
adapted  (cf.  Figure  4.3). 

4.4.4  Effect  of  Initial  Control  Policy,  r(0) 

Figures  4.10  and  4.11  show  the  importance  of  the  initial  con¬ 
trol  on  the  adaptive  control  response.  Figure  4.10  compares  the  model 
response  with  adaptive  control  responses,  with  open-loop  control 
initially  (i.e.  r(0)  =  0)  and  with  positive  feedback  on  the  first  effect 
level,  W-j ,  initially.  Both  runs  were  for  the  same  disturbance  and 
adaptive  loop  gains.  The  response  is  much  more  satisfactory  with  no 
control  initially.  It  is  encouraging  that  the  MRAC  technique  can  recover 
from  a  disasterous  initial  control  policy  such  as  one  involving  positive 
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FIGURE  4.9  SIMULATED  EFFECT  OF  ADAPTIVE  LOOP  GAINS  WITH  FEW 
ELEMENTS  OF  r  ADAPTING 
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FIGURE  4.10  SIMULATED  EFFECT  OF  THE  INITIAL  CONTROL  POLICY  ON 
MRAC,  5^.  =  10  (+50%  peed  PLOW-) 
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feedback. 

Figure  4.11  illustrates  a  logical  initial  control  policy,  well- 
tuned  proportional  gains  for  multiloop  control  (i.e.  several  single 
variable  loops  [5,  Chapter  4])  with  W-j  controlled  by  B-j ,  W2  controlled 
by  B2>  and  C2  controlled  by  S.  The  Appendix  presents  the  multi  loop  con¬ 
trol  matrix  used  for  initial  control  in  Figure  4.11.  Multiloop  control 
as  the  initial  condition  on  r  shows  improved  control  over  initial  open- 
loop  control,  as  would  be  expected.  Figure  4.12  favourably  compares 
adaptive  control  to  multi  loop  proportional  feedback  control  without 
adapting.  Thus,  the  MRAC  technique  can  be  used  to  evolve  a  multivariable 
feedback  control  law  which  improves  on  the  initial  multi loop  control 
policy.  Furthermore,  MRAC  provides  a  systematic  method  for  the  transi¬ 
tion  (during  normal  operation)  from  a  multi loop  to  a  multivariable 
control  technique.  MRAC  also  provides  a  promising  approach  for  tuning 
multi  loop  or  multivariable  control  constants. 

4.4.5  Effect  of  Unseen  Disturbances 

Figure  4.13  demonstrates  that  control  deteriorates  only 
slightly  if  there  is  a  process  disturbance  that  is  not  observed  by  the 
model.  This  robust  feature  of  the  MRAC  technique  is  desirable  since, 
in  practical  situations,  not  all  process  disturbances  are  monitored. 

4.4.6  MRAC  with  a  Physically  Impossible  Reference  Model 

Figure  4.14  shows  a  simulation  run  that  used  a  reference  model 
which  maintained  all  model  states  at  their  initial  value  despite  dis¬ 
turbances  (i.e.  "perfect  control",  x^t)  =  0).  The  gain  of  ^  =  10 
displays  oscillatory  results,  but  a  gain  of  1  exhibits  a  smooth,  although 
sluggish,  response.  This  indicates  that  the  larger  the  model  error,  the 
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FIGURE  4.1.1  SIMULATED  MULTILOOP  CONTROL  AS  AN  INITIAL  CONTROL 
POLICY  =  10  C+50%  FEED  FLOW) 
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FIGURE  4.12 


SIMULATED  COMPARISON  OF  MULTILOOP  AND  MULTIVARIABLE 
MRAC  (+50%  FEED  FLOW) 
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FIGURE  4.13  SIMULATED  EFFECT  OF  UNSEEN  DISTURBANCE  ON  MRAC 
(+20%  FEED  FLOW,  +10%  FEED  COMPOSITION) 
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FIGURE  4.14  SIMULATED  EFFECT  OF  PHYSICALLY  IMPOSSIBLE  REFERENCE 

MODEL,  x  =  0 
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lower  the  adaptive  loop  gain  must  be  to  prevent  oscillatory  responses. 

The  robustness  of  the  adaptive  control  scheme  is  again  demonstrated 
since  even  a  physically  impossible  reference  model  resulted  in  a  less 
than  disasterous  control  scheme. 

4.4.7  Adaptive  Control  in  the  Presence  of  a  Continual  Disturbance 

Figures  4.15  and  4.16  indicate  that  control  gradually  improves 
when  a  continual  disturbance  repeatedly  upsets  the  process.  Figure  4.16 
is  a  continuation  of  Figure  4.15  after  approximately  75  hours  on  MRAC 
and  demonstrates  almost  perfect  model  following.  This  greatly  increases 
the  usefulness  of  adaptive  control  since  control  will  continue  to  improve 
as  disturbances  continue  to  upset  the  process.  These  simulation  results 
support  the  theoretical  conclusion  that  the  error  between  process 
response  and  model  response  eventually  goes  to  zero. 

4.4.8  State  Weighting 

With  all  the  £..  elements  identical,  the  adaptive  control  law, 

'  J 

y- •  =  x  .  e^  P.  5. .  (Chapter  Three),  minimizes  the  errors  between  the 

*  J  r  J  •  *  J 

model  state  variables  and  the  process  state  variables  equally.  However, 
in  practical  situations,  some  state  variables  are  less  critical  than 
others  and  errors  between  these  process  and  model  states  can  be  more 
easily  tolerated.  When  the  optimal  feedback  control  for  the  reference 
model  was  chosen,  the  more  important  states  (W-j ,  W2>  and  C^)  were  weighted 
in  the  quadratic  performance  index  to  obtain  the  optimal  control  law  [5]. 
Although  this  partially  provided  for  state  weighting,  another  approach 
to  state  weighting  was  attempted  by  selecting  ■  values  on  an  individual 

'  J 

basis  to  stress  the  important  states.  Figure  4.17  demonstrates  the  effect 
of  weighting  C2  (i.e.  Xp^)  more  heavily  by  increasing  the  values  of  the 
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FIGURE  4.15  SIMULATED  RESULTS  OF  THE  EFFECT  OF  A  CONTINUING  DISTURBANCE  ON  MRAC  E. . 

t  =  0  TO  150  MINUTES  (CONTROLS  B1 ,  B2,  S  NOT  MONITORED)  ’J 
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FIGURE  4.16  SIMULATED  EFFECT  OF  A  CONTINUOUS  DISTURBANCE  ON 
MRAC ,  t  + 
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FIGURE  4.17  SIMULATED  EFFECT  OF  STATE  WEIGHTING  WITH  ADAPTIVE 
LOOP  GAINS 
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adaptive  gains  in  the  fifth  row  of  s  while  keeping  the  other  .  ele- 

—  '  J 

ments  constant.  It  should  be  noted  that  the  control  of  C ^  improves 
while  control  of  the  other  states  stays  the  same  or  deteriorates  slightly. 
This  state  weighting  method  was  also  demonstrated  in  Figure  4.9  when  all 
the  elements  did  not  adapt.  In  the  simulated  run  in  Figure  4.9,  the 
adaptive  loop  gains  in  the  rows  corresponding  to  C-j  and  h^  were  set 
equal  to  zero  which  gave  essentially  no  weighting  to  these  states.  The 
control  of  these  states  deteriorated  more  notably;  however,  the  improve¬ 
ment  of  the  more  heavily  weighted  states  due  to  the  high  gains  used  in 
their  loops  can  also  be  seen. 

4.4.9  Final  Process  Response  and  Final  Feedback  Matrix 

The  theory  for  MRAC  presented  in  Chapter  Three  stated  that  the 
model  and  process  responses  eventually  become  equal  and  the  process 
feedback  control  matrix  becomes  constant  for  all  values  of  adaptive  loop 
gain.  However,  the  theory  does  not  require  the  final  closed-loop  pro¬ 
cess  matrix,  Ap»  to  be  equal  to  the  closed-loop  model  matrix,  A^,  because 
the  Liapunov  function  total  time  derivative,  V,  is  only  negative- 
semidefinite  rather  than  negative-definite.  Figures  4.15  and  4.16  demon¬ 
strate  that  the  process  and  model  responses  eventually  become  equal.  The 
first  two  matrices  in  Table  4.1  compare  the  model  feedback  matrix  and  the 
final  process  feedback  matrix  (for  adaptive  loop  gains  of  1  and  10)  and 
demonstrate  that  the  eventual  closed-loop  process  matrix  is  not  equal 
to  tne  model  matrix  although  the  error  between  model  and  process  responses 
reduces  to  zero. 

When  an  adaptive  loop  gain  of  100  was  used  with  the  same 
<  iitial  conditions  and  repeated  disturbance  (i.e.  open-loop  initial 
conditions  and  a  40%  step  change  in  feedflow  every  45  minutes),  the  third 
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TABLE  4.1 

MODEL  AND  PROCESS  FEEDBACK  MATRICES 

1)  Optimal  Model  Feedback  Matrix: 


10.78 

-1.607 

-4.818 

0.0 

-19.57 

^FBm 

5.352 

.3598 

.5471 

0.0 

12.49 

7.515 

1.273 

.1833 

24.61 

32.69 

2)  Final ^  Process  Feedback  Matrix  (5. .  =  1  or  10,  Open-Loop  Initial 

■  J 

Control ) 


4.274 

-1.863 

.942 

.625 

-1.782 

3.831 

-  .4540 

.919 

-.018 

.225 

7.477 

-  .6290 

-.1703 

2.217 

1.394 

o 

3)  Approximate  Final  Process  Feedback  Matrix  (5..  =  100,  Open-Loop 
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Initial  Control) 


200,000 

-6000 

150 

30,000 

-9000 

^FBp  = 

160,000 

-3000 

150 

30,000 

-4000 

_  400,000 

1000 

150 

90,000 

-8000__ 

4) 

Final ^  Process 

Matrix  (5..  = 

■  J 

100,  Multi  loop  Initial 

Control , 

Small 

3 

Disturbances  ) 

7.214 

-6.274 

-2.656 

2.118 

-14.60" 

=  FBp  = 

7.543 

-.0192 

-1.564 

2.182 

.6802 

14.44 

2.933 

-5.183 

25.52 

10.90_ 

1. 

After  100  hours  adapting  (elements  almost 

stationary) . 

2. 

Average  of  two 

values  taken 

20  hours  apart  after  200 

hours  adapting. 

3.  4%  disturbance  in  feed  flow  every  20  minutes. 
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feedback  matrix  in  Table  4.1  was  obtained  after  200  hours  of  adapting. 

The  large  gains  in  this  matrix  produced  oscillatory,  stable  control 
with  a  response  similar  to  Figure  4.6  and  the  control  variables 
exhibited  "bang-bang"  (i.e.  open-shut)  behavior.  Therefore,  the  simu¬ 
lation  results  indicate  that  for  certain  conditions,  the  adaptive  con¬ 
trol  law  yields  a  nonstationary  control  matrix  with  large  elements  and 
the  error  between  process  and  model  states  never  reduces  to  zero.  It 
was  also  discovered  that  the  magnitude  of  the  elements  of  the  final 
process  feedback  matrix  could  be  reduced  and  the  ultimate  control 
improved  if  the  initial  control  policy  was  improved,  the  magnitude  of  the 
disturbances  was  reduced,  or  the  adaptive  loop  gains  were  lowered.  For 
example,  a  multi  loop  initial  control  policy  and  small  disturbances  (4% 
step  change  in  feedflow  every  20  minutes  or  40%  changes  every  1000 
minutes)  produced  the  fourth  matrix  presented  in  Table  4.1.  When  the 
initial  control  was  an  optimal  multivariable  policy  based  on  an  inaccurate 
model  (the  first  and  second  effect  holdups  used  in  the  model  were  in 
error  by  50%  and  20%),  an  adaptive  loop  gain  of  1000  produced  a  stationary 
Kpg  with  elements  similar  to  the  optimal  model  feedback  matrix  and  a 
reduction  of  the  error  to  zero. 

Several  differences  exist  between  the  theory  and  the  applica¬ 
tion  of  MRAC  that  could  explain  the  large  final  feedback  matrix  obtained 
in  some  of  the  runs.  The  pseudo-inverse  used  to  generate  the  process 
feedback  matrix,  Kpg,  from  the  adaptive  control  matrix,  r,  (needed  in  all 
processes  where  the  number  of  control  variables  is  less  than  the  number  of 
state  variables)  provides  only  an  approximate  solution  since  the  15 
elements  of  KpB  are  determined  from  the  25  elements  of  r.  This  is  one 
reason  why  there  are  several  possible  final  control  matrices  depending  on 
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the  initial  policy,  adaptive  loop  gain,  and  type  of  disturbance.  The 
pseudo-inverse  approximation  could  be  eliminated  if  a  third-order  refer¬ 
ence  model  was  used.  A  third-order  model  is  available  [5]  and  the 
future  investigation  of  MRAC  using  this  model  should  be  considered  in 
order  to  eliminate  the  pseudo-inverse  problem. 

A  second  difference  between  theory  and  application  is  that  the 
theory  does  not  make  allowance  for  constrained  control  variables; 
however,  the  hybrid  simulation  restricted  the  controls  to  vary  only 
about  the  values  t  1.0  in  perturbated  variables.  Combinations  of  poor 
initial  control  policies,  violent  distrubances ,  and/or  large  adaptive 
loop  gains  violated  these  constraints  and  could  cause  the  control  matrices 
to  become  large  as  shown  by  the  results  presented  in  Table  4.1.  Since  all 
practical  processes  have  constrained  control  variables,  this  is  an 
important  limitation  of  the  theory. 

Finally,  the  theory  was  developed  for  a  continuous  process  and 
reference  model;  however,  the  application  of  the  theory  used  a  discrete 
model  and  sampled  process  states.  The  backward  difference  approximation 
used  for  the  adaptive  algorithm  (i.e.  Equation  (3.48))  could  result  in 
large  r  matrices  (and  hence  a  large  Kpg)  due  to  numerical  instability 
problems.  A  large  time  increment  (sampling  interval),  poor  initial  con¬ 
ditions  (initial  control  policy),  or  large  adaptive  loop  gains  could 
cause  the  numerical  stability  problems.  To  test  this  hypothesis,  the 
sampling  interval  of  64  seconds  was  doubled  and  then  redoubled  with  an 
adaptive  loop  gain  of  10,  an  open-loop  initial  control  policy,  and  40% 
step  changes  in  feed  flow  every  45  minutes.  When  the  sampling  interval 
(and  model  calculations)  were  doubled,  a  constant  KpB  consisting  of 
small  elements  resulted.  However,  when  the  sampling  interval  was  doubled 
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again,  the  process  feedback  matrix  blew  up.  This  supports  the  conclusion 
that  numerical  problems  were  involved  and  is  a  potential  problem  in  the 
discrete  application  of  this  MRAC  technique.  A  higher  order  difference 
approximation  for  Equation  (3.24)  or  a  shorter  sampling  (adapting) 
period  would  correct  this  problem  and  allow  higher  loop  gains;  however, 
additional  calculation  time  would  be  required. 

4.5  CONCLUSIONS 

The  simulation  results  demonstrate  the  ability  of  model  refer¬ 
ence  adaptive  control  to  improve  upon  the  initial  control  of  a  process. 
The  rate  of  adapting  is  determined  by  the  adaptive  loop  gains,  e.  ..,  in  a 

*  J 

manner  similar  to  that  of  the  controller  gains  used  in  a  conventional 
control  scheme. 

All  elements  of  the  r  matrix  need  not  adapt  to  provide  improved 
control,  but  the  guarantee  of  stability  is  then  lost.  For  the  evaporator 
application,  higher  individual  loop  gains  are  allowable  when  fewer 
elements  adapt.  Also,  care  must  be  taken  when  eliminating  adapting 
loops  or  unstable  control  may  result  since  the  open-loop  evaporator  is 
unstable. 

Adaptive  control  can  recover  from  poor  initial  control  policies 
such  as  positive  feedback.  Furthermore,  the  controller  gains  used  in 
standard  single  variable  control  loops  (i.e.  multiloop  control)  provide 
a  good  starting  point  for  adaptive  multivariable  control.  This  could 
make  adaptive  control  an  attractive  method  of  transferring  a  process  from 
multi  loop  control  to  multivariable  control.  The  reference  model  does  not 
have  to  be  an  exact  dynamic  model  of  the  actual  process  but  only  defines 
the  desired  closed-loop  process  response;  therefore,  modelling  errors  are 
not  as  serious  in  adaptive  control  as  they  can  be  in  other  schemes.  MRAC 
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can  also  be  used  for  the  on-line  tuning  of  a  multivariable  scheme  which 
is  based  on  an  inaccurate  model.  The  results  in  Figures  4.13  and  4.14 
demonstrating  good  control  even  though  model  inaccuracies  were  great 
(unseen  disturbances  and  a  physically  impossible  reference  model) 
support  this  conclusion. 

In  the  evaporator  application  under  certain  conditions,  large 
nonstationary  feedback  control  matrices  resulted  and  these  large  matrices 
caused  bang-bang  control  with  the  process  and  model  responses  never 
becoming  equal.  The  conditions  that  caused  the  unstable  results  were 
large  adaptive  loop  gains,  poor  initial  control  policies,  and/or  large 
disturbances.  The  discrete  solution  of  the  adaptive  control  differential 
equation  was  one  reason  for  the  large  KpB  matrix.  Other  differences 
between  the  theory  and  the  application  of  MRAC,  the  pseudo-inverse, 
approximation  and  the  violation  of  control  constraints,  could  also  result 
in  large  feedback  matrices. 

State  weighting  can  be  accomplished  by  using  large  elements  in 
the  rows  of  the  adaptive  loop  gain  matrix  s  corresponding  to  critical 
state  variables.  By  varying  the  individual  values  of  s. .,  it  may  be 
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possible  to  use  higher  adaptive  gains  in  certain  loops  than  are  possible 
when  all  loop  gains  are  equal. 

An  important  result  is  that  control  continues  to  adapt  with 
every  disturbance,  even  those  which  the  reference  model  is  not  aware 
of.  This  would  indicate  that  if  a  sufficiently  low  loop  gain  was  chosen, 
adaptive  control  should  be  continued  after  it  has  improved  the  initial 
control,  since  later  disturbances  will  tend  to  continually  adjust  and 
improve  the  control. 
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CHAPTER  FIVE 

MODEL  REFERENCE  ADAPTIVE  CONTROL  APPLIED  TO 
A  DOUBLE  EFFECT  EVAPORATOR:  EXPERIMENTAL  RESULTS 

5.1  INTRODUCTION 

This  chapter  presents  an  experimental  application  of  model 
reference  adaptive  control  (MRAC)  based  on  Liapunov's  direct  method. 

The  experimental  apparatus  consisted  of  a  pilot  plant  double  effect 
evaporator  that  is  interfaced  to  an  IBM  1800  data  acquisition  and  con¬ 
trol  computer  in  the  Department  of  Chemical  and  Petroleum  Engineering 
of  the  University  of  Alberta.  The  theoretical  background  for  this 
control  method  and  simulation  results  for  the  evaporator  have  been 
presented  in  the  previous  two  chapters.  This  chapter  includes  a  des¬ 
cription  of  the  experimental  control  system  and  of  the  programming 
techniques  used  on  the  IBM  1800.  The  experimental  results  are  then 
presented  and  compared  to  the  simulation  results.  Finally,  conclusions 
concerning  the  feasibility  of  model  reference  adaptive  control  for 
industrial  applications  are  offered. 

5.2  CONTROL  SYSTEM 

The  steps  required  to  implement  model  reference  adaptive 
control  are  presented  in  block  diagram  form  in  Figure  5.1.  An  IBM  1800 
data  acquisition  and  control  computer  is  interfaced  to  a  pilot  plant 
double  effect  evaporator  through  analog  to  digital  converters,  digital 
to  analog  converters,  and  required  transducers.  Figure  5.2  presents  a 
schematic  diagram  of  the  evaporator.  Further  details  concerning  the 
evaporator  equipment  and  instrumentation  can  be  found  in  the  theses  of 
Newell  [3],  Fehr  [4]  and  Andre  [5]. 
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FIGURE  5.1  SCHEMATIC  DIAGRAM  OF  MODEL  REFERENCE  ADAPTIVE  CONTROL 
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FIGURE  5.2  PILOT  PLANT  DOUBLE  EFFECT  EVAPORATOR  WITH  A  CONVENTIONAL  MULTILOOP  CONTROL  SYSTEM 
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The  IBM  1800  incorporates  a  time-shared  version  of  an  industrial 
Direct  Digital  Control  (DDC)  package  [6]  that  can  implement  conventional 
multi  loop  control  on  the  evaporator  through  six  single  control  loops  and 
four  cascaded  control  loops.  However,  in  multivariable  model  reference 
adaptive  control,  the  DDC  control  package  is  used  only  to  access  the 
state  variables  (through  measurements  stored  in  the  DDC  loops)  and  to 
implement  control  (through  setpoint  manipulation  in  the  DDC  loops). 

A  standard  multivariable  control  program  developed  by  Newell 
[3]  executed  the  various  steps  shown  in  Figure  5.1.  Since  one  of  the 
states,  C-j ,  is  not  measured,  the  control  program  performed  the  necessary 
state  estimation  and  also  provided  filtering  of  the  DDC  measurements. 

It  also  performed  the  units  conversion  necessary  to  change  the  measured 
state  variables  into  normalized  perturbation  form  prior  to  the  control 
calculations,  and  to  change  the  calculated  control  signals  from  perturba¬ 
tion  form  to  engineering  units  suitable  for  transmission  to  the  process. 
This  program  was  loaded  into  core  from  disk  and  executed  at  a  standard 
interval  of  64  seconds  in  response  to  a  high  priority  timer  interrupt. 

The  results  from  the  previous  control  interval  were  needed  for  the  cal¬ 
culations  and  were  also  obtained  from  disk.  Because  of  the  slow  disk 
operations  (0.5  seconds  average  access  time),  the  program  required  from 
four  to  seven  seconds  to  execute.  If  the  control  program  had  been  core 
resident,  the  time  required  for  control  calculations  would  have  been  an 
order  of  magnitude  lower.  All  the  control  calculations,  state  estimation, 
filtering  and  units  conversion  were  done  in  FORTRAN  with  only  simple 
Assembler  programs  required  for  communications  between  the  multivariable 
program  and  DDC  (via  the  Process  Variable  Table).  A  complete  listing  of 
the  standard  multivariable  control  program  developed  by  Newell  is 
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available  in  another  research  report  [3]. 

The  computer  program  just  described  is  modular  in  design  and 
only  the  input  and  control  sections  of  the  standard  multivariable  pro¬ 
gram  were  changed  to  implement  multivariable  model  reference  adaptive 
control.  The  subroutines  that  calculated  the  model  states  and  the 
adaptive  control  were  identical  to  those  used  in  the  hybrid  simulation 
studies.  Thus,  the  time  required  to  de-bug  the  real-time  program  was 
greatly  reduced  by  using  a  hybrid  computer  in  the  simulation  studies. 
More  details  of  the  programs,  listings,  logic  flowsheets  and  execution 
for  the  adaptive  control  coreloads  can  be  found  in  the  adaptive  control 
user's  manual  [7]. 

5.3  EXPERIMENTAL  RESULTS 

Fifteen  experimental  runs  were  made  on  the  double  effect 
evaporator  to  investigate  the  MRAC  design  parameters  and  to  confirm  the 
conclusions  drawn  from  the  simulation  results  of  Chapter  Four.  The 
experimental  runs  investigated  the  effect  of  adaptive  loop  gain, 
initial  control  policy,  state  weighting,  repeated  disturbances,  number 
of  adapting  elements,  violation  of  control  constraints,  and  the  opera¬ 
tion  of  the  evaporator  under  MRAC  at  process  conditions  considerably 
different  from  those  used  to  obtain  the  linear  model. 

The  following  figures  show  the  response  of  the  actual  evapora¬ 
tor  to  feed  disturbances  of  approximately  20%,  and  have  an  open-loop 
initial  control  policy  unless  otherwise  specified.  Figure  5.3  illus¬ 
trates  a  typical  evaporator  experimental  run.  Three  of  the  measured 
states  (state  C-j  is  not  measured)  are  presented  on  the  left  half  of  the 
plot  with  T-j ,  first  effect  temperature,  representing  the  first  effect 
enthalpy,  h-|(i.e.  Xpg),  and  the  control  variables  and  observable 
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disturbances  given  in  the  right  half. 

5.3.1  Improvement  in  Control 

Figure  5.3  demonstrates  that  in  spite  of  the  unfavourable 
initial  control  policy,  MRAC  provides  fairly  rapid  model -following  and 
adequate  control.  Figure  5.4  compares  another  experimental  run  with 
different  adaptive  loop  gains  to  the  optimal  feedback  reference  model, 
and  demonstrates  the  improved  control  after  a  second  20%  disturbance 
in  feed  flow. 

Figure  5.5  compares  the  experimental  run  illustrated  in 
Figure  5.3  with  the  simulation  results  for  the  same  disturbance.  As 
illustrated  in  this  figure,  the  experimental  results  are  in  general 
agreement  with  the  simulated  results  except  for  W^,  the  response  of 
this  variable  being  too  oscillatory  initially.  This  could  be  due  to 
process  noise  as  compared  to  the  virtual  absence  of  noise  in  the  hybrid 
simulation.  However,  evaporator  noise  levels  are  quite  small  [8]. 
Modelling  errors  and  the  nonlinear  nature  of  the  evaporator  are  more 
likely  the  cause  of  these  differences.  Figure  5.6  also  compared  simu¬ 
lated  and  experimental  results  and  a  generalization  seems  possible;  W-j 
shows  improved  experimental  control  while  the  control  of  C2  and  W2  is 
better  in  the  simulation.  This  is  also  demonstrated  in  Figure  5.7 
which  compares  experimental  and  simulated  errors  in  W-j .  The  experimental 
error  reduces  to  zero  more  rapidly.  Figure  5.8  demonstrates  that  the 
adapting  parameter  y^ ,  the  element  of  r  that  is  affected  most  by  W-| , 
adapts  more  quickly  in  the  experimental  run.  The  modelling  errors 
between  the  actual  nonlinear  evaporator  and  the  linear  open-loop  model 
used  in  the  simulation  probably  caused  these  differences. 
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FIGURE  5.3.  EXPERIMENTAL  RESULTS  OF  MRAC,  ADAPTIVE  LOOP  GAIN 
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FIGURE  5.4  EXPERIMENTAL  MRAC  AND  REFERENCE  MODEL  (£..  =  0.25, 

c  =  10,  £..  =  1  for  all  other  i,  j) 
i  J 


LEVEL  IN  IN-  H0D  CONCN  IN  PC  GLYCOL 


123 


■ 


CONCN  IN  PC  GLYCOL 


124 


TIME  IN  MINUTES 

FIGURE  5.7  EXPERIMENTAL  W,  ERROR  VS.  SIMULATED  W,  ERROR  (?.,  =  . 
5i5  =  10,  all  ither  Si j  =  1 )  1 
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FIGURE  5.8  EXPERIMENTAL  y51  ELEMENT  VS.  SIMULATED  y51  ELEMENT 
(?i3  =  0.25,  ^i5  =  10,  all  other  z-.  =  1) 
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5.3.2  Effect  of  Adaptive  Loop  Gain 

Figure  5.9  illustrates  adaptive  control  when  =  10.  This 
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plot  demonstrates  more  rapid  adapting  to  the  reference  model  than  for 
the  previous  loop  gains,  however,  oscillations  result.  Figure  5.10  is 
a  continuation  of  Figure  5.9  and  demonstrates  that  continued  oscillations 
occur  after  a  second  step  change  in  flow.  The  state  variables  have  a 
satisfactory  response  with  the  oscillations  being  of  low  magnitude; 
however,  the  control  variables  (B-| ,  B2  and  S)  exhibit  sustained  oscilla¬ 
tions  with  large  amplitudes.  In  practical  situations,  these  fluctuations 
in  control  variables  would  be  undesirable  due  to  excessive  wear  on  equip¬ 
ment  such  as  control  valves.  This  would  indicate  that  £••  =  10  is  above 

*  J 

the  upper  limit  of  the  adaptive  loop  gain  for  satisfactory  control  of  the 
evaporator.  The  feedback  control  matrix  resulting  from  this  run  was  not 
as  large  as  the  optimal  feedback  matrix  used  in  the  reference  model  cal¬ 
culations;  therefore,  the  gain  of  10  did  not  cause  the  numerical  problems 
associated  with  the  large  gains  in  the  simulation.  By  contrast,  the 
simulation  runs  showed  that  =  10  gave  satisfactory  control.  This 
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result  confirms  a  conclusion  derived  from  the  simulation  studies  (Chapter 
Four):  there  is  an  upper  limit  on  the  adaptive  loop  gains  for  the 
evaporator  and  gains  above  this  limit  give  an  oscillatory  response. 

The  disturbances  in  feed  temperature,  most  noticeable  in 
Figures  5.9  and  5.10  and  present  in  all  experimental  runs,  are  a  result 
of  changes  in  the  steam  flow  (S)  to  the  first  effect.  An  interaction 
results  because  the  same  manifold  supplies  steam  to  control  the  process 
and  to  pre-heat  the  feed. 

5.3.3  State  Weighting 


Figure  5.11  demonstrates  the  effect  of  state  weighting  in 
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FIGURE  5.9  EXPERIMENTAL  MRAC 
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FIGURE  5.10  EXPERIMENTAL  MRAC :  SECOND  DISTURBANCE  FOR 
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FIGURE  5.11  EXPERIMENTAL  EFFECTS  OF  STATE  WEIGHTING  AND  REPEATED  DISTURBANCES  (£.-  =  0.25, 
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adaptive  control  and  the  improvement  in  control  for  repeated  disturbances. 
Figure  5.11  used  the  same  adaptive  loop  gains  and  was  subjected  to  the 
same  disturbance  as  in  Figure  5.4  and  demonstrates  the  reproducibility 
of  the  experimental  results.  The  adaptive  loop  gains  were  chosen  as 
£.|3  =  0.25,  £.j5  =  10  and  all  other  =  1  to  weight  the  errors  in  the 
second  effect  concentration,  C2  (i.e.  xp5),more  strongly  than  the  first 
effect  enthalpy,  h-j  (i.e.  xp3).  The  adapting  law  (Chapter  Three) 


Yij  '  xpj  -  -i  5ij  (l’  j  1 . 5) 

can  be  weighted  towards  the  state  by  increasing  the  jt'1  column  of 
the  adaptive  loop  gain  matrix  5.  Simulation  studies  in  Chapter  Four 
demonstrate  that  improved  C2  control  results  with  this  weighting  method 
and  a  comparison  of  the  initial  response  in  Figure  5.3  with  the  initial 
response  in  Figure  5.11  for  the  states  C2  and  T-j  (i.e.  h-j)  demonstrates 
improved  control  in  C2  and  a  deterioration  of  control  in  T-j . 


5.3.4  Effect  of  Repeated  Disturbances 

Figure  5.11  also  demonstrates  that  control  and  adapting  to  the 
reference  model  will  improve  with  continuing  process  disturbances.  The 
initial  response  to  the  first  feed  disturbance  is  poor  because  of  the 
initial  open-loop  control  on  the  system.  When  the  second  step  change  in 
feed  flow  is  made,  control  is  much  better.  Figure  5.7  demonstrates  the 
improved  reference  model  following  for  the  first  effect  level,  W-j . 

These  runs  confirm  the  theoretical  and  simulation  results  that  demonstrated 
continuing  adapting  and  reduction  of  the  error  between  the  model  states 
and  the  process  states  to  zero  with  repeated  disturbances.  Since  the 
steam  control  causes  continuing  upsets  in  the  feed  enthalpy  due  to  the 
common  steam  manifold,  this  result  is  important  for  this  particular  process. 
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5.3.5  Initial  Control  Policy 

Figure  5.12  illustrates  a  good  initial  control  policy,  r(0), 
for  MRAC,  namely,  well -tuned  multi  loop  control  values  for  the  control 
loops  illustrated  in  Figure  5.2.  The  initial  multiloop  control  values 
are  the  same  ones  used  in  the  simulation  runs  with  multiloop  initial 
control  and  are  presented  in  the  Appendix  for  Chapter  Four.  The  initial 
response  is  much  better  than  the  corresponding  run  with  initial  open- 
loop  control  shown  in  Figure  5.11.  The  response  does  not  show  as  much 
improvement  over  the  multiloop  control  as  the  simulation  results  demon¬ 
strated.  However,  the  control  continued  to  adapt  with  MRAC  and  the 
second  disturbance  has  less  effect.  Conventional  multi  loop  control 
(without  adapting)  would  give  essentially  the  same  type  of  response  for 
the  second  disturbance  as  it  did  for  the  first  disturbance. 

5.3.6  Number  of  Adapting  Elements  of  r 

All  elements  of  r  need  not  adapt  for  adequate  control  to  be 
maintained.  In  the  run  depicted  in  Figure  5.13,  only  11  elements  of 
the  r  elements  adapted,  the  same  elements  used  in  the  simulation  results 
of  Figure  4.9  (Chapter  Four).  These  loops  were  chosen  such  that  adapting 
would  be  maintained  for  the  largest  elements  in  the  optimal  KpB  matrix 
used  in  the  reference  model.  The  control  here  is  exceptional  because 
the  loop  gain  chosen,  =  10,  is  probably  close  to  the  best  adaptive 
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gain  for  this  system.  (Note  that  £.  •  =  10  produced  oscillatory  results 

'  J 

when  all  the  gains  adapted.  Figure  5.9.)  This  phenomenon  was  also  demon¬ 
strated  in  the  simulation  results. 

5.3.7  Effect  of  Unseen  Disturbances 


The  effect  of  unseen  disturbances  is  important  in  the  real 
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FIGURE  5.12  EXPERIMENTAL  MRAC  WITH  MULTILOOP  INITIAL  CONTROL  VS.  MULTILOOP  CONTROL  =  0.25 
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process  since  changes  in  feed  composition,  Cp^),  cannot  be  detected 
during  an  experimental  run.  Figure  5.14  illustrates  the  small  effect  a 
disturbance  in  feed  composition,  Cp,  of  10%  not  observed  by  the  model 
has  on  control.  This  was  also  demonstrated  in  the  simulation  results. 

5.3.8  Violation  of  Control  Constraints 

The  experimental  runs  encountered  problems  not  experienced 
with  the  simulation  studies:  the  evaporator  is  nonlinear  (as  compared 
to  the  linear  open-loop  model  used  in  the  simulation  runs),  one  of  the 
state  variables  (C-j )  is  not  measured,  and  there  is  a  physical  constraint 
on  one  of  the  controls  that  is  more  severe  than  the  constraint  used  in 
the  simulation  studies  (i.e.  ±1.0  in  perturbated  control  variables). 

In  this  latter  instance,  the  steam  control  variable,  S,  was  not  able  to 
achieve  a  maximum  flow  rate  of  4.0  Ib/min  that  was  assumed  in  the  simula¬ 
tion  study.  In  fact,  a  flow  rate  of  only  3.0  lb/min  could  be  achieved. 

In  Figure  5.15,  the  +20%  change  in  feed  flow  caused  the  steam  constraint 
and  the  constraint  on  the  second  effect  bottoms  flow  (negative  flow  was 
required  by  the  control  calculations)  to  be  violated  and  in  turn, 
caused  the  control  to  deteriorate  initially.  However,  after  the  control 
has  adapted  sufficiently,  the  control  constraints  are  not  violated  and  a 
much  improved  experimental  response  results. 

5.3.9  Change  in  Operating  Conditions 

Finally,  the  process  was  tested  at  operating  conditions  far 
from  where  the  linearized  model  was  obtained.  A  feed  rate  of  3.5  pounds 
per  minute  was  used  instead  of  the  normal  5  pounds  per  minute.  This 
caused  a  higher  concentration  in  the  second  effect  product  than  is 
normally  produced  since  the  steam  flow  rate  was  kept  constant  at  the 
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FIGURE  5.14  EXPERIMENTAL  MRAC :  EFFECT  OF  UNSEEN  DISTURBANCE  (£..  =  1,  -10%  FEED  COMPOSITION 
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nominal  value  of  2.0  pounds  per  minute.  Figure  5.16  demonstrates  that 
MRAC  improves  control  greatly  and  produced  a  new  feedback  matrix  to 
control  the  process  at  an  unusual  set  of  steady  state  conditions. 

5.4  CONCLUSIONS 

The  experimental  results  supported  all  of  the  conclusions 
drawn  from  the  simulation  studies.  MRAC  performed  satisfactorily  for 
a  wide  range  of  initial  control  policies,  adaptive  loop  gains,  steady 
state  operating  conditions  and  modelling  errors  (such  as  unseen  distur¬ 
bances).  The  value  of  the  adaptive  loop  gains,  is  important  in 
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determining  the  type  of  response  obtained.  Low  values  of  5..  produced 
more  sluggish  responses  while  larger  values  resulted  in  faster,  more 
oscillatory  responses.  Thus  adaptive  loop  gains  and  conventional  con¬ 
troller  gains  exhibit  analogous  behavior. 

All  runs  demonstrated  continuously  improving  control  when 
repeated  disturbances  acted  on  the  process.  This  would  indicate  that 
although  a  low  value  of  £. .  can  produce  initially  sluggish  control, 

*  J 

this  control  will  continue  to  improve  and  provide  closer  model -fol lowing 
with  repeated  disturbances. 

The  large  feedback  matrices  produced  in  the  simulation  results 
were  not  encountered  in  the  experimental  runs.  However,  the  experimental 
runs  were  not  conducted  for  the  same  length  of  time  as  the  simulation 
runs  nor  was  an  adaptive  loop  gain  above  10  used.  However,  with  longer 
runs  and  higher  loop  gains,  large  control  matrices  would  probably  occur 
in  the  experimental  runs,  also. 

A  well -tuned  multi  loop  control  system  provided  an  excellent 
initial  control  policy.  The  adaptive  control  scheme  also  provided 
adequate  control  for  unseen  disturbances  which  is  an  important  evaporator 
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FIGURE  5.16  EXPERIMENTAL  MRAC  AT  LOW  FEED  RATES 
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consideration  since  changes  in  feed  composition  are  undetected.  Also, 
all  elements  of  r  need  not  adapt  for  the  process  response  to  follow  the 
model  response  in  a  satisfactory  fashion. 

The  experimental  results  demonstrated  that  adaptive  control 
can  survive  constraints  on  the  control  variables  that  are  not  included 
in  the  model  calculations  if  the  adaptive  loop  gain  is  small  enough. 

The  unsuspected  steam  constraint  caused  poorer  initial  control  when  a 
disturbance  caused  the  violation  of  this  constraint;  however,  later 
disturbances  tended  not  to  violate  the  constraint  as  control  improved. 

Adaptive  control  demonstrated  further  insensitivity  to  modell¬ 
ing  errors  when  it  maintained  adequate  control  at  a  steady  state  far 
from  that  at  which  the  process  was  modelled. 

Two  practical  applications  of  model  reference  adaptive  control 
designed  by  Liapunov's  direct  method  may  be  postulated.  This  technique 
could  prove  useful  in  tuning  multivariable  control  systems  since  it  does 
not  require  an  accurate  process  model.  For  example,  optimal  control 
schemes  depend  greatly  on  model  accuracy  for  their  success.  If  the 
fixed  parameter  control  scheme  derived  from  a  model  fails  to  perform 
satisfactorily,  MRAC  could  tune  the  multivariable  constants  instead  of 
requiring  new  model  calculations.  Also,  MRAC  provides  a  systematic 
approach  for  transferring  a  process  from  a  standard  multi  loop  control 
configuration  to  multivariable  control. 
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CHAPTER  SIX 
CONCLUSIONS 

The  purpose  of  this  work  was  to  investigate  a  model  reference 
adaptive  control  (MRAC)  system  designed  by  Liapunov's  direct  method  of 
stability  analysis.  The  theory  for  proportional  feedback,  feedforward 
model  reference  adaptive  control  as  developed  by  Winsor  and  Roy  [1]  and 
Porter  and  Tatnall  [2]  was  extended  to  include  integral  and  setpoint 
control  modes.  Proportional  feedback  MRAC  was  then  applied  to  a  linear 
model  of  a  double  effect  evaporator.  The  work  included  44  simulation 
runs  on  an  EAI  590  hybrid  computer  in  order  to  test  the  various  design 
parameters  of  the  MRAC  approach.  The  simulation  runs  investigated  the 
effects  of  adaptive  loop  gain,  initial  control  policy,  unmeasured 
disturbances,  state  weighting,  repeated  disturbances,  and  the  number 
of  adapting  elements.  Fifteen  experimental  runs  were  then  made  on  the 
evaporator/IBM  1800  control  configuration  to  confirm  the  simulation 
results.  Since  the  actual  process  was  nonlinear  and  had  an  unexpected 
constraint  on  steam  flow  rate,  the  experimental  runs  also  tested  addi¬ 
tional  aspects  of  the  control  scheme.  Table  6.1  presents  a  summary  of 
the  results  and  conclusions  of  the  previous  three  chapters.  The  follow¬ 
ing  two  sections  present  recommendations  for  the  industrial  application 
of  MRAC  and  suggestions  for  future  work. 

6.1  INDUSTRIAL  APPLICATIONS  OF  MODEL  REFERENCE  ADAPTIVE  CONTROL 

Many  types  of  processes  could  benefit  from  a  MRAC  system. 
Highly  nonlinear  processes  that  are  subject  to  changing  operating  condi¬ 
tions  can  be  difficult  to  control  with  a  fixed  parameter  control  scheme. 
Modern  multivariable  control  schemes  are  usually  based  on  a  linear 


■ 

' 


SUMMARY  OF  DESIGN  CONSIDERATIONS  IN  MRAC  BY  LIAPUNOV'S  DIRECT  METHOD 


142 


143 


-o 

cu 


o 

a 


vo 


co 

<c 


144 


process  model.  MRAC  could  be  used  to  tune  the  fixed  parameter  control 
scheme  if  process  model  inaccuracy  or  changing  process  conditions  caused 
unsatisfactory  results. 

For  application  of  model  reference  adaptive  control  to 
industrial  processes,  it  is  recommended  that  small  adaptive  loop  gains 
and  a  good  initial  control  policy  be  used.  The  MRAC  system  designed  by 
Liapunov's  direct  method  for  linear  systems  is  theoretically  stable; 
however,  process  nonlinearities  and  numerical  problems  that  arise  with 
the  discrete  form  of  the  adaptive  control  law  can  cause  stability 
problems.  To  eliminate  the  numerical  problems,  smaller  sampling 
(adapting)  intervals  or  higher-order  difference  equations  could  be  used 
to  solve  the  adaptive  control  equations  (3.24,  3.25),  but  a  sacrifice 
in  calculation  time  and  computer  storage  would  result. 

The  control  matrices  should  be  monitored  if  a  poor  initial 
control  policy  is  used  in  order  that  a  combination  of  severe  disturbances 
or  changing  operating  conditions  does  not  result  in  large  control 
matrices  and  "bang-bang"  control.  Intermittent  application  of  MRAC,  as 
opposed  to  continuous  application,  could  be  used  on  an  industrial  pro¬ 
cess  in  order  to  save  computer  time.  After  MRAC  has  improved  the 
initial  control  and  the  process  and  model  states  have  become  almost 
equal,  the  MRAC  algorithm  calculations  could  be  suspended  until  a  set- 
point  change  is  made  that  alters  the  process  operating  conditions.  To 
demonstrate  this  application,  MRAC  was  terminated  on  the  simulated  pro¬ 
cess  after  forty-five  minutes  and  the  resulting  constant  control  matrix 
was  used  to  control  the  process  thereafter  for  a  similar  disturbance. 
Figure  6.1  demonstrates  the  performance  of  the  fixed  parameter  control 
system  that  resulted  with  5..  =  10  and  open-loop  initial  control. 
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The  reference  model  could  be  more  easily  derived  than  the 
optimal  closed-loop  model  used  in  this  study.  Less  complex  design 
techniques  that  do  not  rely  on  dynamic  programming  calculations  could 
be  used  to  specify  the  model.  Matrix  assignment  [3]  or  eigenvalue 
assignment  [4,5],  where  the  desired  type  of  closed-loop  response  is 
specified,  would  be  ideal  for  model  reference  adaptive  control  since 
the  closed-loop  model  matrix  is  selected  by  the  design  engineer  and 
does  not  require  complex  calculations.  In  this  study,  adequate  control 
was  achieved  even  when  all  the  model  state  variables  remained  unchanged 
despite  disturbances  (i.e.  the  assignment  of  infinitely  negative 
eigenvalues)  in  the  simulation  studies  (Chapter  Four).  This  result 
indicates  that  MRAC  is  feasible  with  any  arbitrarily  assigned, 
asymptotically  stable,  reference  model. 

6.2  FUTURE  WORK 

The  simulation  and  experimental  results  have  indicated  a 
number  of  areas  that  need  further  investigation  in  the  design  and 
application  of  MRAC  systems. 

6.2.1  Feedforward,  Integral,  and  Setpoint  Adaptive  Control 

The  addition  of  a  feedforward  control  configuration  can  greatly 
improve  the  response  of  many  processes  compared  to  feedback  control 
alone.  Furthermore,  the  variations  in  the  disturbance  coefficient 
matrix,  D,  due  to  changing  disturbance  vector  steady-states,  are  often 
greater  than  variations  in  the  process  matrix.  A,  and  adaptive  feed¬ 
forward  control  could  partially  compensate  for  these  changes.  Similarly, 
since  integral  action  is  usually  desirable  to  eliminate  offset,  the 
application  of  the  theory  developed  in  this  work  for  integral  adaptive 
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control  could  find  extensive  use.  Setpoint  control  is  essential  in 
industrial  processes  and  because  setpoint  changes  vary  process  operating 
conditions  and  hence  the  linear  process  model,  a  further  investigation 
of  setpoint  adaptive  control  should  be  conducted. 

6.2.2  Selection  of  Adaptive  Loop  Gains 

The  adaptive  loop  gains  were  found  to  behave  analogously  to 
standard  controller  gains  and  had  an  appreciable  effect  on  the  adaptive 
control  response.  An  a  priori  method  of  determining  a  good  value  of 
these  gains  would  be  desirable.  Gromyko  and  Sankovskii  [6]  have  also 
demonstrated  the  effect  of  adaptive  loop  gain  in  simulation  studies 
and  have  developed  a  complicated  method  for  determining  the  optimal 
value  of  adaptive  loop  gain  from  the  model  and  process  matrices  if  the 
variations  in  process  parameter  are  known.  Further  investigation  in 
the  area  of  adaptive  loop  gains  is  needed. 

6.2.3  Selection  of  Reference  Model 

Many  possible  asymptotically  stable  reference  models  could  be 
chosen  for  a  specific  process.  Because  a  reasonably  accurate  open- 
loop  model  was  available  for  the  evaporator,  optimal  feedback  control 
was  used  in  selecting  the  closed-loop  reference  model  for  this  study. 
Other  reference  models  can  be  chosen  on  the  basis  of  eigenvalue  assign¬ 
ment  [5]  and  these  reference  models  should  also  be  tested.  The 
evaporator  is  a  nonlinear  process  and  control  of  the  evaporator  using  a 
fifth-order  model  essentially  tested  MRAC  with  a  different  order  pro¬ 
cess  than  reference  model.  Loan  [7,8]  has  recently  investigated  process 
identification  by  Liapunov's  direct  method  for  both  lower-order  and 
higher-order  models.  The  theory  he  has  developed  could  be  investigated 
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on  the  evaporator  by  using  reference  models  other  than  the  fifth-order 
model  used  in  this  study  and  determining  the  effect  of  reference  model 
order  on  control . 

6.2.4  Constrained  State  and  Control  Variables 

The  MRAC  theory  outlined  in  this  work  does  not  consider  con¬ 
straints  on  the  state  or  control  variables.  Gromyko  and  Sankovskii 
[6]  have  used  constrained  control  variables  in  their  paper  on  adaptive 
loop  gains.  They  imply  that  stability  is  guaranteed  only  for  a  certain 
range  of  adaptive  loop  gains  if  there  are  constraints  on  the  variables. 
The  further  implications  of  constrained  variables  should  be  investigated 
since  all  processes  operate  within  certain  limits. 

6.2.5  Selection  of  the  P  Matrix  Used  in  the  Adaptive  Control  Law 

The  rows  of  the  Liapunov  matrix  P  appear  in  the  adaptive 
control  law 

T 

v..  =  x  •  e  P .  £ . . 

PJ  -  nj 

and  have  a  large  effect  on  the  adaptive  control  response.  The  P  matrix 
is  obtained  from  the  Liapunov  equation 

AmT  P  +  P  Am  =  -  Q 

=m  =  =  =m  2 

and  its  value  depends  on  an  arbitrarily  specified  positive-definite 
matrix  Q  and  the  closed-loop  model  matrix  Am-  Therefore,  another  design 
option  is  available  when  the  Q  matrix  is  specified.  The  effect  of  this 
design  parameter  on  control  and  its  possible  uses  (e.g.  state  weighting) 
should  be  studied. 
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6.2.6  Discrete  Approximation  of  the  Adaptive  Control  Law 

The  numerical  stability  problems  encountered  in  the  evaporator 
application  could  be  eliminated  with  higher-order  difference  approxima¬ 
tions  to  solve  the  adaptive  control  algorithms.  For  example,  the 
feasibility  of  doing  a  separate  analysis  on  the  adaptive  algorithm 

T 

y . .  =  x  .  e  P.  £ .  . 

YiJ  PJ  -  -i 

to  determine  the  maximum  adaptive  loop  gains  for  a  specific  P  matrix, 
difference  equation,  and  control  interval  should  be  considered. 

6.3  CONCLUSIONS 

Model  reference  adaptive  control  designed  from  Liapunov 
stability  considerations  has  been  applied  to  a  double  effect  evaporator. 
In  both  simulation  and  experimental  studies,  MRAC  proved  useful  for 
changing  from  multiloop  control  to  multivariable  control.  The  scheme 
was  able  to  recover  from  very  poor  initial  control  strategies.  It  was 
used  to  tune  a  multivariable  control  system  based  on  an  inaccurate 
model.  Therefore,  this  MRAC  scheme  could  be  used  in  industrial  situa¬ 
tions  where  it  is  desired  to  change  from  multiloop  to  multivariable 
control  or  to  tune  existing  multivariable  control  schemes.  Feedback 
and  setpoint  control  modes  should  be  adequate  for  most  multivariable 
MRAC  systems;  however,  feedforward  and  integral  adaptive  control  could 
be  added  to  processes  where  offsets  are  not  desirable. 

The  implementation  of  model  reference  adaptive  control  is 
straightforward,  however,  the  nonlinear  differential  nature  of  the 
control  law  and  the  extensive  matrix  operations  can  cause  calculation 
time  and  storage  requirements  to  be  an  order  of  magnitude  larger  than 
for  a  similar  fixed  parameter  control  scheme.  If  calculation  time 
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and/or  storage  requirements  are  limited,  adapting  only  certain  critical 
loops  or  the  intermittent  application  of  model  reference  adaptive  con¬ 
trol  (e.g.  with  setpoint  changes)  could  be  considered. 


NOMENCLATURE  FOR  CHAPTER  TWO 


Alphabetic 

A  heat  transfer  area 

=  state  coefficient  matrix 

B  control  coefficient  matrix 

C  concentration 

Cp  heat  capacity 

e  error 

E  activation  energy 

F  feed  flow 

h  level  or  local  heat  transfer  coefficient 

H  heat  of  reaction 

J  performance  criterion 

k  frequency  factor 

K  controller  gain 

L  load  or  length 

m  controller  output 

n  concentration 

P  perimeter 

q  flow 

r  rate  of  reaction 

R  setpoint,  resistance  to  flow,  or  gas  constant 

s  Laplace  variable 

t  time 

T  temperature,  sampling  period 

U  overall  heat  transfer  coefficient 
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NOMENCLATURE  FOR  CHAPTER  TWO  (continued) 


control  vector 

v  volume  or  river  velocity 

V  vol ume 

x  composition  or  spacial  coordinate 

_x  state  vector 

z  Z-transform  variable 

a  constant  of  diffusion 

3  time  scaling  factor 

A  incremental  quantity 

e  error 

p  density 

Tj  integral  time 

t  time  delay  or  time  constant 

w  frequency 

Subscripts 

c  coolant  or  controller 

v  val  ve 

Superscipts 


time  derivative 
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NOMENCLATURE  FOR  CHAPTER  THREE 

Alphabetic 

state  coefficient  matrix 
positive  constant 
control  coefficient  matrix 
first  effect  bottoms  flow 
second  effect  bottoms  flow 
first  effect  concentration 
second  effect  concentration 
feed  concentration 
output  coefficient  matrix 
disturbance  vector 
disturbance  coefficient  matrix 
error  vector 
feed  flow 

controller  parameters  in  M.I.T.  MRAC  technique 
first  effect  enthalpy 
feed  enthalpy 
identity  matrix 
control  matrix 
state  dimension 
disturbance  dimension 
positive -definite  symmetric  matrix 
positive -definite  symmetric  matrix 
r  control  dimension 

S 


At 


steam  flow 

discrete  time  interval 
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NOMENCLATURE  FOR  CHAPTER  THREE  (continued) 


u 

V 

W1 

W, 

x 

y 

z 


control  vector 
Liapunov  function 
first  effect  holdup 
second  effect  holdup 
state  vector 

output  vector  of  integrated  states 
"integral"  state  vector 


Greek 


A 


ij 


ij 


Xij 


V 

K 

t 

ft 

0) 


IJ 


’IJ 


IJ 


state  difference  matrix,  Am  -  Ap 
elements  of  A 

adaptive  feedback  control  matrix 
elements  of  r 

discrete  control  coefficient  matrix 

discrete  disturbance  coefficient  matrix 

disturbance  difference  matrix,  Dm  - 

=m  =p 

elements  of  A 

adaptive  feedforward  loop  gains 
adaptive  feedback  loop  gains 
discrete  state  coefficient  matrix 
adaptive  feedforward  control  matrix 
elements  of  ft 


Subscripts 


discrete 


FB 


feedback 
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NOMENCLATURE  FOR  CHAPTER  THREE  (continued) 


FF  feedforward 

I  integral 

m  model 

p  process 

sp  setpoint 

ss  steady  state 

Superscripts 

k  control  interval 

T  transpose 

*  perturbation  variables  or  augmented  vector 

time  derivative 
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NOMENCLATURE  FOR  CHAPTER  FOUR 


Alphabetic 

A  closed-loop  state  matrix 

B1  first  effect  bottoms  flow 

B2  second  effect  bottoms  flow 

C-j  first  effect  concentration 

C2  second  effect  concentration 

Cp  feed  concentration 

d  disturbance  vector 

e  error  vector 

F  feed  flow 

h-j  first  effect  enthalpy 

KpB  feedback  control  matrix 

P  positive-definite  matrix 

S  steam  flow 

u  control  vector 

V  Liapunov  function 

W-j  first  effect  holdup 

W2  second  effect  holdup 

x  state  vector 

Greek 

a.j j  elements  of  A=  ^  -  Ap 

r  adapting  feedback  matrix 

y •  •  elements  of  r 

'  J  = 

h  adaptive  loop  gain  matrix 

adaptive  loop  gains,  elements  of 
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NOMENCLATURE  FOR  CHAPTER  FOUR  (continued) 


Subscripts 

m  model 

p  process 

Superscripts 


T 


transpose 
time  derivative 
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NOMENCLATURE  FOR  CHAPTER  FIVE 


Alphabetic 

B-j  first  effect  bottoms  flow,  lbs/min 

B2  second  effect  bottoms  flow,  lbs/min 

C.|  first  effect  concentration,  weight  percent 

C2  second  effect  concentration,  weight  percent 

d  disturbance  vector 

e  error  vector 

F  feed  flow,  lbs/min 

hp  feed  enthalpy,  BTU/min 

h-j  first  effect  enthalpy,  BTU/min 

0-j  first  effect  overhead,  lbs/min 

02  second  effect  overhead,  lbs/min 

P-j  first  effect  pressure,  psia 

P2  second  effect  pressure,  psia 

P.  column  of  Liapunov  Matrix  P 

J  “ 

SIFB  steam  flow,  lbs/min 

S  steam  flow 

T-j  first  effect  temperature 

Tp  feed  temperature,  °F 

u  control  vector 

x  state  vector 

Greek 

r  adapting  control  matrix 

E  adaptive  loop  gain  matrix 
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NOMENCLATURE  FOR  CHAPTER  FIVE  (continued) 


adaptive  loop  gains 

Superscripts 

T  transpose 

*  time  derivative 


Subscripts 


P 


probess 


NOMENCLATURE  FOR  CHAPTER  SIX 


Alphabetic 

A  state  equation  coefficient  matrix 

D  state  equation  coefficient  matrix 

e  error  vector,  model  state  -  process  state 

P  Liapunov  matrix,  positive-definite 

Q  arbitrary  positive-definite  matrix 

x  state  vector 

Greek 

.  adaptive  loop  gains 

*  J 

adapting  control  elements 

Superscripts 

time  derivative 
T  transpose 

Subscripts 

p  process 

m  model 
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APPENDIX  FOR  CHAPTER  THREE 
NUMERICAL  VALUES  OF  MATRICES  USED  FOR  MRAC 

a)  Optimal  Continuous  Feedback  Matrix,  KpB,  for  the  Control  Law  u(t) 

=  =FB  SpC*) 

The  optimal  KpB  was  obtained  from  the  IBM  360-67  program  MARIC^ 

with  the  state  weighting  matrix  Q  =  diagonal  [10,1,1,10,100],  the  control 

r  T  T 

weighing  matrix  R  =  0,  in  the  performance  index,  J  =  (x  Q  x  +  u  Ru)  dt, 


25.87 

-3.254 

-11.18 

-.1143 

-54.39 

-FB 

16.57 

.4828 

1.200 

-12.48 

74.99 

7.218 

1.169 

.0201 

29.05 

30.95 

b)  Continuous  Model  Matrices 

The  ^pB  given  above  was  substituted  into  xm  = 

(A  +  B 

+  (or  ^  H  £ + 

B  KpB)  to  yield 

“  -1.269 

-.0381 

-.2174 

.9558 

-5.744 
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-.0755 

.1255 

0.0 

0.0 

Sm  ■ 
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-.0074 

-.0502 

-2.098 

4.780 

-.6854 
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1.  DACS  Centre  Program  Library. 
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c)  Discrete  Model  Matrices 

The  open-loop  evaporator  model,  xp  =  A  xp  +  B  u  +  D  d,  can  be 
written  in  discrete  form  as,  xp(k+l)  =  $  xp(k)  +  A  u(k)  +  e  d(k),  where 


t  = 

[  efi*  dt 

J  0 

A  = 

[T  e&V-Q 
Jo 

B  dt 

0  = 

fT 

J  0 

D  dt 

For  T  =  64  seconds,  the  matrices  are 

"1.0 

.0008 
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.9223 
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0 

1  = 
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.0042 

.4396 
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.0001 

_0 

.0391 
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0 

A  = 

.1569 

0 
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-.0138 

.0848 

-.0406 

_  .0137  - 

.0432 

0 

"  .1182 

0 

-.0050” 

-.0351 

.0785 

.0049 

0  = 

-.0136  - 

.0002 

.0662 

.0012 
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-.0058 

L -  .0019 

.0016 

.0058 

The  open-loop  discrete  evaporator  model 


xp(k+l)  =  $  xp(k)  +  A  u(k)  +  0  d(k) 
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with  feedback  control 


y(k)  =  =FBd  *p(k) 

yields  the  closed-loop  reference  model  used  in  the  MRAC  scheme. 

xm(k+l)  =  C|  +  A  KpBd)  xm(k)  +  e  d(k) 

xm( k+1 )  =  <|>  x  (k)  +  0  d(k) 

-m  =m  -m  =m  - 


The  optimal  discrete  feedback  matrix,  KpBd,  obtained  from  GEMSCOPE  [22] 
with  the  same  weighting  matrices  as  in  the  continuous  case  and  T  =  64 
seconds  is: 


“10.78 

-1.607 

-4.818 

0.0 

-19.57" 

=  FBd 

5.352 

.3598 

.5471 

0.0 

12.49 

_  7.515 

1.273 

.1833 

24.61 

32.69_ 

and  the  closed-loop  matrices  that  result  are 
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.0662 

=m 

.0012 

0 

-.0058 

-.0019 

.0016 

.0058 

. 
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d)  Solution2  of  aJ”  P  +  P  Am  =  -  g  for  P 


3.224 

1.070 

-.0483 

.1258 

-5.444 

1.070 

5.367 

.0430 

.0038 

-1 .968 

- 

-.0483 

.0430 

.1629 

-.0103 

.0352 

.1258 

.0038 

-.0103 

.3529 

.2327 

-5.444 

-1.968 

.0352 

.2327 

10.34 

e)  Results  Using  the  Pseudo-Inverse  of  r 

Equation  (3.4)  gave  r(t)  =  B  KpB C t )  but  r  is  solved  from  a  back¬ 
ward  difference  equation  using  the  discrete  matrix  A 


(A)-1 


r(k)  =  a  kpBd(k) 
kFBdCk)  *  CaV1  aT  r(k) 


-.3303 

.4618 

6.258 

0 

.6241 

= 

-9.545 

-.0206 

.  -.2785 

0 

-5.087 

_-19 .80 

-.1992 

-2.699 

-24.61 

-10.82 

2.  Solved  by  the  method  described  by  Chen,  C.F.  and  Sheih,  L.S.,  "A  Note 
on  Expanding  P  A  +  A^P  =  -Q",  IEEE  Trans,  on  Auto.  Control,  AC-13,  4, 
448  (1968). 


■  . 
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APPENDIX  FOR  CHAPTER  FOUR 
ADAPTING  ELEMENTS  FOR  RUNS  WITH  A  REDUCED  NUMBER 
OF  ADAPTING  ELEMENTS 


FIGURE 

NUMBER 


4.7 


4.7 


4.7 


4.8 


4.9 


NUMBER  OF  ELEMENTS 

ADAPTING  ADAPTING  ELEMENTS  OF  r 

21  All  elements  except 

y34,y54 

11  First  and  fifth  rows  of  r  and 

y44 

7  Yn  »y44,y15,y25,y35,y45,y55 

4  Y11 *y44*y15*y55 

7  Yn  ,y44  ,Y1 5  *y25  *y35  ,y45  *y55 


FEEDBACK  CONTROL  MATRIX  FOR  MULTILOOP  CONTROL 


0  0  0  0  -4.89 

3.52  OOO  0 

000  15.8  0 


■ 

■ 


. 


• 
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PROCESS  STEADY  STATE  OPERATING  CONDITIONS 

x  -  State  Vector 

W-|  -  First  effect  holdup  -  45.5  lb. 

C-j  -  First  effect  concentration  -  4.59  %  glycol 

h-|  -  First  effect  solution  enthalpy-189.2  BTU/lb. 

W2  -  Second  effect  holdup  -  41.5  lb. 

C2  -  Second  effect  concentration  -  10.11  %  glycol 

-  Control  Vector 

S  -  Steam  flowrate  to  the  first  effect  -  2.0  Ib/min 
B-j-  First  effect  bottoms  flow  -  3.48  lb/mi n 

B2~  Second  effect  bottoms  flow  -  1.58  Ib/min 

d  -  Disturbance  Vector 

F  -  Feed  flow  -  5.0  lb/mi n 

Cp  -  Feed  concentration-3.2  l  glycol 
hp  -  Feed  enthalpy  -  156.9  BTU/lb. 


' 

' 


